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ABSTRACT 

We have compiled a pseudo-snapshot data set of two-colour observations from the Hubble 
Space Telescope archive for a sample of 53 rich LMC clusters with ages 10^ — 10^" yr. We 
present surface brightness profiles for the entire sample, and derive structural parameters for 
each cluster, including core radii, and luminosity and mass estimates. Because we expect the 
results presented here to form the basis for several further projects, we describe in detail the 
data reduction and surface brightness profile construction processes, and compare our results 
with those of previous ground-based studies. 

The surface brightness profiles show a large amount of detail, including irregularities in 
the profiles of young clusters (such as bumps, dips, and sharp shoulders), and evidence for 
both double clusters and post core-collapse (PCC) clusters. In particular we find power-law 
profiles in the inner regions of several candidate PCC clusters, with slopes of approximately 
—0.7, but showing considerable variation. We estimate that 20 ± 7 per cent of the old cluster 
population of the LMC has entered PCC evolution, a similar fraction to that for the Galactic 
globular cluster system. In addition, we examine the profile of R136 in detail and show that it 
is probably not a PCC cluster. 

We also observe a trend in core radius with age that has been discovered and discussed in 
several previous publications by different authors. Our diagram has better resolution however, 
and appears to show a bifurcation at several hundred Myr. We argue that this observed rela- 
tionship reflects true physical evolution in LMC clusters, with some experiencing small scale 
core expansion due to mass loss, and others large scale expansion due to some unidentified 
characteristic or physical process. 
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1 INTRODUCTION 

The star cluster system of the Large Magellanic Cloud (LMC) 
is unique in containing rich star clusters of masses compara- 
ble to Galactic globular clusters, but covering a wide age range 
(10*^ - 10^" yr) and being close enough for detailed observation. 
It therefore offers a seemingly perfect opportunity for studies of 
all aspects of star cluster astronomy, from cluster formation, evo- 
lution and dynamics, to luminosity and mass function studies, as 
well as investigations which shed light on the evolution history of 
the entire LMC. It is surprising then, that while catalogues of high 
resolution surface brightness profiles and structural parameters ex- 
ist for Galactic globular clusters (e.g., Trager, King & Djorgovski 
(1995)), no such uniform catalogue exists for a large sample of the 
rich LMC clusters. 

There has however been some considerable activity in this 
field. In particular, significant numbers of surface brightness and/or 
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density profiles have been published for young and intermediate 
age clusters (Elson, Fall & Freeman (1987), hereafter EFF87; El- 
son, Freeman & Lauer (1989); Elson (1991), hereafter E91; El- 
son (1992)); old clusters (Mateo 1987); and clusters with different 
spatial distributions - in the LMC disk (Kontizas, Chrysovergis & 
Kontizas (1987a)); in the LMC disk and within 5 kpc of the rota- 
tion centre (Chrysovergis, Kontizas & Kontizas (1989)); and in the 
LMC halo (Kontizas, Hadjidimitriou & Kontizas (1987b); Metaxa, 
Kontizas & Kontizas (1988)). These studies are all ground-based 
and therefore suffer from problems - primarily crowding and see- 
ing related ~ which limit their resolution, particularly in the inner 
regions of clusters. This in turn renders the derivation of key param- 
eters, such as the core radius, rather uncertain. In addition, each 
set of authors uses a different data set and makes different mea- 
surements (e.g., Elson and collaborators construct l/-band surface 
brightness profiles; Mateo uses B-band profiles; and Kontizas and 
collaborators use number density profiles), so most of the important 
derived values are not strictly comparable between studies. 

We have taken advantage of the presence of a large number of 
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observations of rich LMC clusters in tlie Hubble Space Telescope 
(HST) archive to compile a high resolution data set, unaffected by 
the problems which beset ground-based studies. Rather than simply 
providing an atlas of surface brightness profiles from this data, the 
aim of the present study is to obtain a statistically homogeneous 
set of profiles and key structural parameters for as many clusters 
as possible, for the purposes of differential comparison. To this 
end, we require a sample as free from bias as possible, and we 
have therefore applied a uniform selection and reduction process to 
the available observations. Because these procedures are rather de- 
tailed, and because we expect our measurements to form the basis 
of several additional projects, in this paper we discuss at length the 
data set and its selection (Section 2) and the reduction and profile 
construction processes (Sections 3 and 4 respectively). In Section 
5 we present the surface brightness profiles and structural parame- 
ter measurements, compare these results with those of the authors 
previously mentioned, and examine some of the interesting sub- 
groups present in the sample - such as binary clusters and post 
core-collapse clusters. Finally, in Section 6 we observe and discuss 
the relationship between core radius and age for the LMC cluster 
system. This relationship has previously been studied (Elson et al. 
(1989); E91; Elson (1992)); however our measurements are able to 
provide new insight into this problem. 

The data presented in Tables 1, 2, 4, and 6, and the 
surface brightness profiles (Fig. 6) are available on-line at 
http://www. ast. cam. ac. uk/STELLARP OP S/LMC Musters/. 



2 THE CLUSTER SAMPLE 
2.1 Observations 

The observational basis of this project is the presence in the HST 
archive of a snapshot siffvey of LMC (and SMC) clusters (HST 
project 5475). This data set consists of two Wide Field Planetary 
Camera 2 (WFPC2) exposures per cluster, through the F450W and 
F555W filters respectively, and with exposure times covering the 
range 20-300 s in F555W and 40-600 s in F450W, dependent on 
the age of the cluster under observation. The data were acquired 
between 1994 January 27 and 1994 December 25, but mostly be- 
fore the WFPC2 cool-down on 1994 April 23. 

Upon retrieval, we discovered that the sample was neither as 
large nor as varied as desired, and consequently we returned to the 
archive to obtain observations of as many additional clusters as pos- 
sible. We located a further six suitable studies (HST projects 5114, 
5897, 5904, 5916, 7307, and 8134), consisting of detailed observa- 
tions of several clusters each. To maintain the uniformity of our se- 
lection as far as possible, we attempted to mimic the snapshot data 
by using only two frames per cluster, with exposure times as close 
as possible to the ranges described above, even though for many of 
the clusters from these six studies, an abundance of archival frames 
is available. To avoid biases such as that due to observer selection 
of "interesting targets," it was necessary to limit the use of these 
extra frames to match the least observed clusters. This fit with our 
primary aim of investigating the LMC cluster system as a whole. 
For this purpose we required data for as large a sample of clus- 
ters as possible, but to be reduced within a reasonable time frame. 
Therefore, we had to balance ultimate accuracy (maximum data) 
for any given cluster in the sample against overall reduction time. 
Had our aim been simply to provide an atlas of surface brightness 
profiles, or to study one or two clusters in great detail, the extra data 
would certainly have been included. Nonetheless, we were careful 



not to unnecessarily degrade the quality of our data. In cases where 
the neglect of additional data would have compromised our results 
(e.g., in extremely crowded clusters) we included minimal numbers 
of extra frames in order to proceed. 

In the end, complete reproduction of the original snapshot data 
was not possible, with the only uniform choice of filters between 
the extra six studies being F555W and F814W. For nine clusters we 
retrieved both long (~ 500 s) and short (~ 10 s) exposures in order 
to overcome data degradation due to severe crowding and/or satura- 
tion (as mentioned above). The final cluster sample is listed in Table 
1 together with the observational details^. The sample consists of 
53 LMC clusters spanning the full age range. The entire sample is 
observed through the F555W filter, providing a unique high reso- 
lution data set. Thirty-one of the clusters are also observed through 
the F450W filter, with the remaining 22 clusters having F814W as 
their second filter - we note however that for our surface brightness 
profile construction process, the choice of the second colour is of 
Uttle consequence. 

2.2 Literature data 

As a supplement to the above data set, and as a reference point 
for this and future publications, we have compiled nomenclature, 
position, age and metallicity data from the literature for the com- 
plete cluster sample. This information is displayed in Table 2. As 
with the data set construction, the emphasis is on obtaining as ho- 
mogeneous a compilation as possible, while still maintaining the 
integrity and accuracy of the data. This means using the results of 
large surveys or studies for the most part, and supplementing these 
with the results of high quality observations where necessary. In 
the cases where several such high quality studies are available, we 
have chosen the values which represent approximately the conver- 
gence of the data. For several clusters, scarcity of information has 
necessitated the selection of older or lower quality data. This com- 
pilation does not purport to be a complete survey of the available 
literature; rather it is intended to provide a consistent set of age and 
metallicity estimates for use in this and subsequent projects. 

2.2.1 Cluster names and positions 

We have taken position and nomenclature information from 
the Simbad Astronomical Database {http://simbad.u-strasbg.fr/). 
While almost all recognized identifiers for a cluster are included in 
this database, we have not compiled all such identifiers for our sam- 
ple, rather only the most common. In general the principal designa- 
tion is an NGC number, with the other identifiers and correspond- 
ing catalogues as follows: Hodge, (Hodge 1960); SL, (Shapley & 
Lindsay 1963); LW, (Lynga & Westerlund 1963). In the Simbad 
database, J2000.0 positions for LMC objects are in general from 
the catalogue of Bica et al. (1999). Using these positions we cal- 
culate the projected angular distance Ropt to the optical centre of 
the LMC at a = 05''20"56^ 5 = -69°28'4l" (J2000.0) (Bicaet 
al. 1996), and the projected angular distance Rrot to the H I rota- 
tion curve centre at a = 05''26"'39^ 6 = -69° 15' 36" (J2000.0) 
(Rohlfs et al. 1984; Westerlund 1990). We later derive more accu- 
rate centres for each cluster; however, the literature values listed 
here are never substantially different from oui calculated positions. 

t Notes to Table 1: 

" The archived images of this cluster are incorrectly labelled as NGC 2156. 
The archived images of this cluster are incorrectly labelled as SL 633. 
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Table 1. Cluster list and observation details. 



Cluster Program Principal Frame Secondary Frame 



Name 


ID 


Filter 


Data set 


Date 


Time (s) 


Filter 


Data set 


Date 


Time (s) 


NGC1466 


5897 


F555W 


u2xj0103t 


22/10/1995 


260 


F814W 


u2xj0101t 


22/10/1995 


260 


NGC1651 


5475 


F555W 


u26m0r02t 


01/02/1994 


120 


F450W 


u26m0r01t 


01/02/1994 


230 


NGC1711 


5904 


F555W 


u2y80502r 


18/10/1997 


300 


F814W 


u2y80504r 


18/10/1997 


300 


NGC1718 


5475 


F555W 


u26m0j02t 


06/02/1994 


120 


F450W 


u26m0j01t 


06/02/1994 


230 


NGC1754 


5916 


F555W 


u2xq0103t 


21/10/1995 


500 


F814W 


u2xq0109t 


21/10/1995 


600 






F555W 


u2xq0101t 


21/10/1995 


20 


F814W 


u2xq0106t 


21/10/1995 


20 


NGC1777 


5475 


F555W 


u26m0u02t 


23/11/1994 


120 


F450W 


u26m0u01t 


23/11/1994 


230 


NGC1786 


5897 


F555W 


u2xj0203t 


19/10/1995 


260 


F814W 


u2xj0201t 


19/10/1995 


260 


NGC1805 


7307 


F555W 


u4ax0206r 


25/07/1998 


140 


F814W 


u4ax020br 


25/07/1998 


300 


NGC1818 


7307 


F555W 


u4ax3005r 


25/09/1998 


140 


F814W 


u4ax300cr 


25/09/1998 


300 


NGC1831 


5475 


F555W 


u26ml002t 


27/01/1994 


40 


F450W 


u26ml001t 


27/01/1994 


80 


NGC1835 


5916 


F555W 


u2xq0203t 


18/10/1995 


500 


F814W 


u2xq0209t 


18/10/1995 


600 






F555W 


u2xq0201t 


18/10/1995 


20 


F814W 


u2xq0206t 


18/10/1995 


20 


NGC1841 


5897 


F555W 


u2xj0708t 


14/11/1995 


800 


F814W 


u2xj0707t 


14/11/1995 


800 


NGC1847 


5475 


F555W 


u26mlb02p 


10/04/1994 


20 


F450W 


u26mlb01p 


10/04/1994 


40 


NGC1850 


5475 


F555W 


u26mle02t 


15/05/1994 


20 


F450W 


u26mle01t 


15/05/1994 


40 


NGC1856 


5475 


F555W 


u26ml302t 


06/02/1994 


30 


F450W 


u26ml301t 


06/02/1994 


60 


NGC1860 


5475 


F555W 


u26ml202t 


07/02/1994 


40 


F450W 


u26ml201t 


07/02/1994 


80 


NGC1866 


8134 


F555W 


u5ay0601r 


21/08/1999 


350 


F814W 


u5ay0604r 


22/08/1999 


350 


NGC1868 


5475 


F555W 


u26m0z02t 


10/04/1994 


100 


F450W 


u26m0z01t 


10/04/1994 


200 


NGC1898 


5916 


F555W 


u2xq0304t 


10/12/1995 


500 


F814W 


u2xq0309t 


10/12/1995 


600 


NGC1916 


5916 


F555W 


u2xq0403t 


10/12/1995 


500 


F814W 


u2xq0409t 


10/12/1995 


600 






F555W 


u2xq0401t 


10/12/1995 


20 


F814W 


u2xq0406t 


10/12/1995 


20 


NGC1984 


8134 


F555W 


u5ay0901r 


08/07/2000 


350 


F814W 


u5ay0904r 


09/07/2000 


350 






F555W 


u5ay0903r 


09/07/2000 


10 


F814W 


u5ay0906r 


09/07/2000 


10 


NGC2004 


5475 


F555W 


u26mld02t 


21/04/1994 


20 


F450W 


u26mld01t 


21/04/1994 


40 


NGC2005 


5916 


F555W 


u2xq0503t 


19/10/1995 


500 


F814W 


u2xq0509t 


19/10/1995 


600 






F555W 


u2xq0501t 


19/10/1995 


20 


F814W 


u2xq0506t 


19/10/1995 


20 


NGC2011 


8134 


F555W 


u5ay0801r 


19/08/1999 


350 


F814W 


u5ay0804r 


20/08/1999 


350 






F555W 


u5ay0803r 


19/08/1999 


10 


F814W 


u5ay0806r 


20/08/1999 


10 


NGC2019 


5916 


F555W 


u2xq0603t 


18/10/1995 


500 


F814W 


u2xq0609t 


18/10/1995 


600 






F555W 


u2xq0601t 


18/10/1995 


20 


F814W 


u2xq0606t 


18/10/1995 


20 


NGC2031 


5904 


F555W 


u2y80301t 


20/10/1995 


300 


F814W 


u2y80304p 


20/10/1995 


300 


NGC2100 


5475 


F555W 


u26mlc02n 


21/04/1994 


20 


F450W 


u26mlc01t 


21/04/1994 


40 


NGC2121 


5475 


F555W 


u26m0x02t 


02/02/1994 


120 


F450W 


u26m0x01t 


02/02/1994 


230 


NGC2136 


5475 


F555W 


u26ml702t 


11/02/1994 


30 


F450W 


u26ml701t 


11/02/1994 


60 


NGC2153'' 


5475 


F555W 


u26ml502t 


19/04/1994 


30 


F450W 


u26ml501t 


19/04/1994 


60 


NGC2155 


5475 


F555W 


u26m0k02t 


01/02/1994 


120 


F450W 


u26m0k01t 


01/02/1994 


230 


NGC2156 


8134 


F555W 


u5ay0301r 


22/08/1999 


350 


F814W 


u5ay0304r 


22/08/1999 


350 






F555W 


u5ay0303r 


22/08/1999 


10 


F814W 


u5ay0306r 


22/08/1999 


10 


NGC2157 


5475 


F555W 


u26ml902t 


26/03/1994 


30 


F450W 


u26ml901t 


26/03/1994 


60 


NGC2159 


8134 


F555W 


u5ay0401r 


16/10/1999 


350 


F814W 


u5ay0404r 


16/10/1999 


350 


NGC2162 


5475 


F555W 


u26m0t02t 


23/08/1994 


120 


F450W 


u26m0t01t 


23/08/1994 


230 


NGC2164 


8134 


F555W 


uSayOlOlr 


24/09/1999 


350 


F814W 


u5ay0105r 


24/09/1999 


350 


NGC2172 


8134 


F555W 


u5ay0502r 


14/10/1999 


350 


F814W 


u5ay0504r 


14/10/1999 


350 


NGC2173 


5475 


F555W 


u26m0m02t 


02/02/1994 


120 


F450W 


u26m0m01t 


02/02/1994 


230 


NGC2193 


5475 


F555W 


u26m0p02t 


30/01/1994 


120 


F450W 


u26m0p01t 


30/01/1994 


230 


NGC2209 


5475 


F555W 


u26m()w02t 


01/02/1994 


120 


F450W 


u26m0w01t 


01/02/1994 


230 


NGC2210 


5897 


F555W 


u2xj0403t 


13/12/1995 


260 


F814W 


u2xj0401t 


13/12/1995 


260 


NGC2213 


5475 


F555W 


u26m0v02t 


25/12/1994 


120 


F450W 


u26m0v01t 


25/12/1994 


230 


NGC2214 


5475 


F555W 


u26ml802t 


16/02/1994 


30 


F450W 


u26ml801t 


16/02/1994 


60 


NGC2231 


5475 


F555W 


u26m0s02t 


06/02/1994 


120 


F450W 


u26m0s01t 


06/02/1994 


230 


NGC2249 


5475 


F555W 


u26m0y02t 


11/08/1994 


120 


F450W 


u26m0y01t 


11/08/1994 


230 


NGC2257 


5475 


F555W 


u26m0d02t 


05/02/1994 


300 


F450W 


u26m0d01t 


05/02/1994 


600 


SL663'' 


5475 


F555W 


u26m0102t 


01/02/1994 


120 


F450W 


u26m0101t 


01/02/1994 


230 


SL842 


5475 


F555W 


u26m0h02t 


01/02/1994 


120 


F450W 


u26m0h01t 


01/02/1994 


230 


SL855 


5475 


F555W 


u26m0q02t 


01/02/1994 


120 


F450W 


u26m0q01t 


01/02/1994 


230 


H0DGE4 


5475 


F555W 


u26m0o02t 


01/02/1994 


120 


F450W 


u26m0o01t 


01/02/1994 


230 


HODGE 11 


5475 


F555W 


u26m0f02t 


01/02/1994 


300 


F450W 


u26m0f01t 


01/02/1994 


600 


HODGE14 


5475 


F555W 


u26m0n02t 


18/09/1994 


120 


F450W 


u26m0n01t 


18/09/1994 


230 


R136 


5114 


F555W 


u2hk0304t 


25/09/1994 


23 


F814W 


u2hk030rt 


25/09/1994 


40 






F555W 


u2hk0302t 


25/09/1994 


3 


F814W 


u2hk030qt 


25/09/1994 


5 
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2.2.2 Cluster ages 

Age determinations for LMC clusters are widely scattered across 
the literature, and involve a large variety of techniques. The 
most useful homogeneous data set we have located is the colour- 
magnitude diagram (CMD) study and age calibration of Geisler 
et al. (1997), which provides age estimates for 15 clusters in our 
sample. Following this lead we have adopted where possible for 
the remaining clusters age determinations based on high resolution 
photometric studies resulting in CMDs. Pre-eminent amongst these 
are HST studies, and we adopt the results of Rich, Shara & Zurek 
(2001) for 5 intermediate age clusters and the results of Olsen et 
al. (1998) for 5 old clusters. We take the age determinations from 
the groimd-based studies of E91 for 9 clusters and those of Dirsch 
et al. (2000) for a further three clusters. For NGC 1805 and NGC 
1818 we adopt the ages from the thorough discussion and review 
of de Grijs et al. (2002b) and similarly for R136 we take the age 
as discussed by Sirianni et al. (2000). Finally, we fill the remaining 
holes in our sample with the age estimates of Elson & Fall (1988), 
who provide an age calibration based on a large sample of litera- 
ture CMDs. We correct these estimates to a distance modulus of 
18.5 using their prescribed method. Wc could not locate a reliable 
age for NGC 1916, but Olsen et al. (1998) do present a CMD for 
this cluster along with CMDs and age estimates for five old LMC 
clusters. The CMD for NGC 1916, though suffering from serious 
differential reddening, is not evidently different from those for the 
other five clusters, so we average the ages of these clusters to obtain 
an estimate for NGC 1916. 



2.2.3 Cluster metallicities 

The definitive study of LMC cluster metallicities is by Olszewski et 
al. (1991) who derive metallicities for ~ 70 LMC clusters from ob- 
servations of the calcium triplet in the spectra of giants. From this 
study we adopt metallicities for 23 of the clusters in our sample. We 
also adopt one metallicity (NGC 1841) from a subsequent paper in 
the same series (Suntzeff et al. 1992). A further three abundances 
(NGC 1850, NGC 2004, NGC 2100) are taken from the spectro- 
scopic study of Jasniewicz and Thevenin (1994). The authors report 
standard deviations of ±0.03 dex in the averages of several obser- 
vations for each of these clusters; however we prefer their errors per 
single observation of ~ 0.2 dex, which are more consistent with the 
other errors in Table 2. We adopt another two abundances from the 
infra-red spectra of Oliva & Origha (1998), and one (NGC 1866) 
from the VLT spectroscopy of Hill et al.(2000). These complete 
the spectroscopic studies suitable for our data set. For a few clus- 
ters we therefore rely on abvmdances derived from CMD fits and 
corresponding to the ages discussed in Section 2.2.2. Four come 
from the study of Dirsch et al. (2000), and two from the HST pho- 
tometry of Rich et al. (2001). For NGC 1805 and NGC 1818 we 
again take the ranges adopted by de Grijs et al. (2002b) after their 
detailed literature review (see also Johnson et al. (2001)), and sim- 
ilarly for R136 we again adopt the abundance from the discussion 
of Sirianni et al. (2000); see also Hunter et al. (1995). For the re- 
maining 14 clusters we were not able to obtain suitable abundance 
measurements from the literature so for each of these we average 
the literature results for all the clusters in our sample of comparable 
age. We emphasize that these values in particular are simply rather 
crude estimates and that any use of them must account for this. In 
this paper, we adopt these abundances solely for the purpose of esti- 
mating mass to light ratios later in the analysis - calculations which 
in any case are rather insensitive to the adopted metallicity. Natu- 



rally, any subsequent use must be carefully judged on the nature of 
the calculation at hand. 



3 PHOTOMETRY 

As part of the archive retrieval process, all frames were reduced 
according to the standard HST pipeline, using the latest available 
calibrations. This is preferable to using the original data as reduced 
soon after observation because the latest calibrations have a much 
longer baseline for calculation and are therefore more accurate. 
Bearing in mind that our data set covers a span of more than five 
years of HST observations, we also maintain the homogeneity of 
our reduction process by using calibrations from a single epoch. 

For photometric measurements we have found Dolphin's 
HSTphot (Dolphin 2000a) to be the most suitable software pack- 
age. Our large number of observations requires the reduction pro- 
cess to be run automatically where possible, without sacrificing 
data integrity. HSTphot is well suited for running in batch mode 
and is specifically tailored to the reduction of WFPC2 frames, rou- 
tinely accounting for the severely under-sampled PCAVFC point 
spread functions (PSFs) and their variation due to sub-pixel posi- 
tioning. In particular, the multiphot routine has proven to be the 
most useful for our data set. This program first aligns accurately 
(to fractions of a WFC pixel) and then performs simultaneous pho- 
tometry on multiple images of the same field, regardless of the filter 
combination. This makes it perfect for the current data set - with 
two frames per field, each through a different filter. Conventionally, 
this would make the removal of cosmic rays and cross identification 
of objects complicated and unreliable. By solving the two frames 
simultaneously, both these problems are accounted for. 

Before starting multiphot each frame is readied with the image 
preparation utilities included with the HSTphot package. This pro- 
cedure includes the masking of bad Aat?^ as flagged by the STScI 
data quality image which accompanies each observation, a first at- 
tempt at the removal of cosmic rays (using a routine based on the 
IRAF task CRREJ), the removal of hot-pixels not flagged by the data 
quality image (using a cr-clipping algorithm), and the robust deter- 
mination of a background image - used by multiphot for stellar 
detection and measurement. 

Photometric measurements are made using multiphot in PSF 
fitting mode. While aperture photometry mode is also an option, we 
have found that PSF fitting produces main sequences on our CMDs 
that are at least as narrow as those from aperture photometry us- 
ing 2 WFC pixel and 3 PC pixel radii. This is presumably because 
multiphot is designed specifically for measurements involving the 
unique WFPC2 PSFs. The detection algorithm used by multiphot 
is similar in principle to that of DoPHOT, being based around the 
recursive location and subtraction of objects at increasingly lower 
threshold levels. We set a minimum threshold level for detection 
of 3cr above the background. In addition, we enable a feature of 
multiphot which calculates an adjustment to the background image 
before each photometry measurement using the pixels just beyond 
the photometry radius (Dolphin 2000a). This is designed to account 
for rapidly varying backgrounds such as those expected in observa- 
tions of the crowded central regions of globular clusters. 

The object classification, sharpness and x parameters pro- 
duced from multiphot' s PSF fitting are used to keep the photom- 

t Such as known bad pixels, saturated pixels, charge traps, bad columns, 
and the vignetted region between chips 
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Table 2. Literature nomenclature, position, age and metallicity data for the cluster sample. 
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12 
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12 
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05''38"*43* 


-69° 06' 03" 
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14 
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Reference Ust: I. de Grijs et al. (2002b); 2. Dirsch et al. (2000); 3. Elson (1991); 4. Elson & Fall (1988); 5. Geisler et al. (1997); 6. HUl et al. (2000); 7. 
Hunter et al. (1995); 8. Jasniewicz & Thevenin (1994); 9. Johnson et al. (2001); 10. Oliva & OrigUa (1998); 11. Olsen et al. (1998); 12. Olszewski et al. 
(1991); 13. Rich et al. (2001); 14. Sirianni et al. (2000); 15. Suntzeff et al. (1992). 
* Calculated metallicity (or age for NGC 1916), as described in the text. 

"* Relative to the optical centre of the LMC bar, at a = 05''20"'56*, S = -69°28'41" (J2000.0) (Bica et al. 1996). 

''Relative to the Hi rotation centre of the LMC, at a = 05''26"*39*, 5 = -69°15'36" (J2000.0) (Rohlfsetal. 1984; Westerlund 1990) 
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etry clean of non-stellar objects and spurious detections. The ob- 
ject classification parameter is a determination as to which proto- 
type PSF the measured profile of an object conforms to - the op- 
tions being stellar (or a marginally resolved stellar pair), extended 
(galaxy), or single pixel (cosmic ray). The sharpness value for an 
object is a measure of how sharp its profile is in comparison with 
the stellar PSF, being negative if the profile is too broad, zero if the 
object is a perfect fit and positive if the profile is too sharp. Finally, 
the X value simply measures the quality of fit of the PSF to the 
object profile. After some experimentation, we found that select- 
ing objects with sharpness between —0.6 and -1-0.6 and x < 3.5 
provided a mostly complete stellar sample for the purpose of con- 
structing surface brightness profiles. For the construction of clean 
CMDs, these values would be somewhat stricter. 

The final photometric measurements obtained from multiphot 
are from PSFs corrected for geometric distortion (Holtzman et al. 
1995a), the filter-dependent plate scale changes (determined em- 
pirically) (Dolphin 2000a), and the 34th row error (Shaklan, Shar- 
man & Pravdo 1995; Anderson & King 1999). The magnitudes are 
also corrected for charge-transfer efficiency (CTE) effects using the 
longest baseline calibration available - that of Dolphin (2000b). 
Notably, this calibration also provides corrections for "warm" data 
- that is, observations taken before the WFPC2 cool-down of 1994 
April 23. This is especially important for our data set, since 25 
clusters were observed in warm conditions. The stellar magnitudes 
are also corrected for PSF residuals and to an aperture of O'.'S us- 
ing groups of isolated stars selected according to a set of strict 
criteria (Dolphin 2000a). Finally, multiphot uses the zero-points 
from the Dolphin (2000b) calibration. We do not convert from the 
HST instrumental system to the Johnson-Cousins system because 
this transformation is not important for the construction of surface 
brightness profiles, instead adding unnecessary scatter to the pho- 
tometry. Dolphin (2000a) estimates the limiting photometric and 
astrometric accuracy of multiphot to be 0.011 to 0.014 mag and 
0.05 pixels (2.5 x 10~^ arcsec on the PC and 5.4 x 10~^ arcsec 
on the WFCs) respectively. 



4 SURFACE BRIGHTNESS PROFILES 

Constructing a surface brightness profile for a given cluster is sim- 
ple in concept; however in practice, our calculations were compli- 
cated by the high resolution and peculiar chip geometry of WFPC2. 
The procedure is similar to that for ground-based data however, and 
we broadly follow the techniques outlined by Djorgovski (1987). 

4.1 Astrometry and centre determination 

One problem presented by the data was due to the four chip struc- 
ture of WFPC2. In order to maintain the accurate relative spatial 
positioning of stars from chip to chip, we overlaid a uniform coor- 
dinate system, accounting for the geometric distortion from HST's 
optical system, the changes in pixel scale from the PC to WFCs 
and the separations between the chips. This was achieved using the 
IRAF STSDAS task METRIC, which converts the chip and pixel co- 
ordinates for a list of stars to pixel coordinates relative to the WFC2 
chip, making the appropriate geometric corrections, and then con- 
verts these corrected pixel coordinates to (a, 6) (J2000.0) using 
the positional information in the image headers. This header in- 
formation is not always reliable, and can introduce errors into the 
{a, S) so calculated. To avoid this, we took the corrected pixel co- 
ordinates from METRIC for each star and derived angular separa- 



tions using the WFC pixel scale of 0.0996 arcsec pixeP^, rather 
than taking the (a, 5) for each star and calculating angular separa- 
tions from these. All centering calculations, annulus construction, 
and completeness and area corrections (see below) used these cor- 
rected pixel coordinates. The only results which incorporate image 
header positional information are the calculated centre coordinates 
expressed in (a, 5) (see Table 4). 

The next task was to locate accurately the centre of any given 
cluster^. Poor centering will tend to artificially flatten a surface 
brightness profile, and this in turn can lead to large systematic er- 
rors in any structural parameters derived from the profile, and may 
obscure important dynamical information such as the presence of a 
post-collapse core. There are several algorithms for locating a clus- 
ter's centre which rely on the symmetry properties of the cluster, 
one example being the mirror-autocorrelation method described by 
Djorgovski (1987). However, such algorithms were not completely 
suitable for our data. Unlike Galactic globular clusters, many Mag- 
ellanic Cloud clusters do not have smooth profiles, but instead are 
clumpy, irregular and not particularly symmetric, making it diffi- 
cult to apply symmetry-algorithms robustly and imiformly across 
the sample. In addition, and compounded by the high resolution 
and small field of view of WFPC2, several clusters have very low 
surface brightness, adding to the problem. The high resolution in 
particular means that it is not sufficient to simply apply spatial 
smoothing to most clusters, because the magnitude of the smooth- 
ing required degrades the intrinsic accuracy of the centre determi- 
nation, and the process therefore loses its value. 

Instead, we employed a Monte Carlo style method to locate 
our cluster centres. In this procedure, each chip is first split into 
boxes of equal area (~ 100 WFC pixels on a side). For each box, 
the surface brightness is calculated by adding up the fluxes for all 
stars in the box (including completeness corrections - see Section 
4.2.2) and dividing by the area of the box. The search is then nar- 
rowed to the region covered by the box with the highest surface 
brightness, and the eight boxes surrounding it. It does not matter if 
this region falls across more than one CCD - wc simply account 
for the vignetted area (where no stars are detected) in the subse- 
quent calculations. If the box with the highest surface brightness 
falls right on the edge of a CCD, where it is not joined by another 
CCD, then the number of surrounding boxes is less (five if the cen- 
tral box is on an edge; three if it lies in a corner) and the amalga- 
mated region is correspondingly smaller. Within this amalgamated 
region, points are randomly generated and the surface brightness 
calculated (as for the boxes) in a circle of radius r about each point. 
If part of one of these circles falls off the edge of a chip or over a 
vignetted region, this "lost area" is accounted for in the calculation. 
After N tries, the point corresponding to the circle with the highest 
surface brightness is taken to be the cluster centre. 

This routine has the disadvantage that it is very good at finding 
the brightest star in a cluster, and although this does represent a sur- 
face brightness peak, it generally does not correspond to the over- 
all cluster surface brightness peak. To avoid this, we would ideally 
like to exclude the few brightest stars in a cluster from the above 
calculations, and instead use a sub-sample which should closely 
trace the cluster's surface brightness profile (i.e., a sample which 
has its greatest density coincident with the cluster's surface bright- 
ness peak). Main sequence stars are such an ensemble, and so for 
each cluster we impose colour and magnitude limits on the cluster 
CMD to select only these stars for the centre determination. In gen- 

§ Or more correctly, the central surface brightness peak. 
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eral this procedure worked well; however for clusters with severe 
central crowding, it was necessary to relax the limits and include 
brighter stars from the giant branch (or upper main sequence in the 
case of very young clusters). 

While inelegant and computationally expensive, our random 
sampling method is robust and easily adapted for clusters of differ- 
ent concentration and richness simply by variation of r and N. As 
a consistency check, we run the centering algorithm independently 
on both colour frames for a field. For N = 2000 the centering is 
repeatable to a limiting accuracy of approximately ±10 WFC pix- 
els - approximately ±1". As noted above however, the conversion 
from pixel coordinates to (a, 5) may introduce large errors if the 
image header information is inaccurate. A comparison of our cal- 
culated centres (Table 4) with those in Table 2 shows good agree- 
ment, meaning that we have not introduced any large errors. We 
note however that this comparison is not sensitive to small (several 
arcsecond) errors. Notwithstanding this, we expect that our calcu- 
lated centres are at least as accurate as the literature values. 

4.2 Annulus construction and flux corrections 

For each frame (two per cluster), four sets of circular annuli were 
constructed about the cluster centre. Two sets had narrow annulus 
widths of l'.'5 and 2" respectively and were designed for sampling 
the central regions of a cluster. To this end, the 1'.'5 annuli were 
calculated to a radius of ~ 20" and the 2" annuli to ~ 30". The two 
other annulus sets had larger widths of 3" and 4" respectively, and 
were used primarily for sampling the outer (less dense) regions of a 
cluster. These two profiles were therefore extended to the maximimi 
radii possible. Although many LMC clusters have slightly elliptical 
isophotes at extended radii, there is little evidence to suggest that 
the mass distribution in inner regions is in general not spherical 
(cf. E91) - hence we fit circular annuli. In addition, because of the 
small WFPC2 field of view, none of the images in the data set fully 
cover any cluster, making the use of elliptical annuli impractical. 

Because of the high resolution of our images we simply 
counted stars to calculate the surface brightness for an annulus. For 
each set, the surface brightness Hi of the i-th annulus is given by: 

3 = 1 

where hi and Oj are the outer and inner radii of the annulus respec- 
tively, Ns is the number of stars in the annulus, and Fj is the flux 
of the j-th star. The factors Ai and Cj are the area correction for 
an annulus and the completeness correction for a star respectively, 
and must be determined before the annulus is constructed. We out- 
line the meaning of these factors and their methods of calculation 
individually below. 

4.2.1 Area corrections 

The area correction Ai for the i-th annulus is used simply to nor- 
malize the flux in the annulus to that for a full annulus. This is 
necessary because the shape of WFPC2 means that for all clus- 
ters, most annuli are not completely covered by the field of view. 
Since the flux through an aimulus is directly related to the annulus 
area, variations in the fractions of annuli covered cause artificial 
fluctuations in a surface brightness profile, and must therefore be 
accounted for. 

The process of determining the area correction for an aimulus 
is complicated by the WFPC2 chip geometry (including the small 



separations between the CCDs), the centering of the cluster on the 
camera and the roll angle of HST at the time of observation. The 
arbitrary nature of the second two factors from cluster to cluster 
once again led us to resort to an inelegant but robust Monte Carlo 
scheme. For a given cluster we determine corrected pixel coordi- 
nates for the sixteen CCD comer pixels which define the WFPC2 
field of view, and from these we derive the border equations of the 
four CCDs. A large number of points are then randomly generated 
over the full area of each annulus, and using the CCD boundary 
equations we determine which points fall on the WFPC2 camera, 
and which do not. The area correction Ai for the i-th annulus is the 
total number of points generated divided by the number "imaged". 
Again, while computationally inefficient, this method is very robust 
and can easily account for the arbitrary geometry of any particular 
observation, with an accuracy limited only by the number of ran- 
dom points. A small amount of experimentation showed ten thou- 
sand points per annulus to be sufficient. To avoid introducing large 
uncertainties into the data we do not use annuli for which the frac- 
tion covered falls below a third. This limits the maximum radius for 
a surface brightness profile to be ~ 75 — 80". 

An example of the process is shown in Fig. 1 - a pointing for 
NGC 1841. Annuli of width 4" have been drawn about the centre 
to 100", and for every second annulus the random points falling in 
the field of view have been plotted. The accuracy with which the 
process handles the complicated geometry is quite evident - note 
particularly the small separations of the chips due to the vignetted 
regions masked by HSTphot, and the significant errors these can 
cause to an annulus of the right radius which falls tangent. At 100" 
the area correction factor is approximately 4, and in practice the 
process was halted at r = 78". 

4.2.2 Completeness corrections 

Even with HST resolution, crowding and saturation can cause sig- 
nificant numbers of stars to be missed by automated detection 
software such as HSTphot. This missing flux can seriously affect 
a cluster's surface brightness profile and must be accounted for. 
To quantify the correction we used the artificial star routine at- 
tached to multiphot. For each cluster, we set this routine to generate 
~ 3.5 X 10® stars per CCD (i.e., ~ 1.4 x 10® stars per cluster) on 
a CMD, the limits of which are set to be two magnitudes above the 
saturation limit and two magnitudes below the faint limit for an im- 
age, and 0.5 magnitudes redder than the reddest region of the clus- 
ter's CMD and 0.5 magnitudes bluer than the bluest region. The 
stars are spatially distributed according to the flux of the image. 
Each star is placed on the image and solved (one at a time) by mul- 
tiphot, using the same settings as for the real data. Stars are flagged 
if they are not recovered. The output photometry file is exactly sim- 
ilar to the multiphot output for real stars, and wc run this artificial 
photometry through our classification, sharpness and x selection 
criteria, again flagging stars which are removed from the sample. 

To determine the completeness function for a given chip, the 
fake stars are binned according to their x- and j/-position on the 
chip, magnitude and colour. Bin sizes are typically 160 pixels in 
X and y in uncrowded regions (independent of chip), 0.2 magni- 
tudes in brightness and 0.25 magnitudes in colour. A large majority 
of our sample suffers little or no crowding, even in the central re- 
gions of a cluster. This means that the binning resolution described 
above, particularly in the positional sense, is perfectly adequate to 
account for the gradual spatial variations in completeness. This is 
demonstrated below, in the example for NGC 2213, which is a typ- 
ical cluster in the sample. For severely crowded regions however. 
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the completeness may suffer rapid and significant spatial variation, 
and in these cases the resolution of the positional binning was in- 
creased to account for this. In the very worst cases, the inclusion 
of short exposures (see Section 4.2.3 below) usually significantly 
alleviated any crowding problems. 

The completeness c for a given bin is the number of success- 
fully recovered (unflagged) stars divided by the total number of 
stars generated in that bin, and the completeness correction Cj for 
the j-\h star in an annulus the inverse of c for the appropriate bin. To 
avoid large uncertainties, we eliminated any stars with complete- 
ness less than cj = 0.25. Experimentation showed that variation of 
this limit (i.e., 0.05 < ci < 0.5) had a negligible effect on most 
profiles, reflecting the fact that incompleteness was not an issue for 
the large majority of clusters. 

Examples of the completeness functions for two clusters (the 
relatively open cluster NGC 2213 and the extremely compact clus- 
ter NGC 2005) are shown in Fig. 2 together with the PC F555W 
images of the centres of these clusters. For ease of visualization 
the completeness functions have been integrated over position and 
colour to be functions of brightness only. For NGC 2213, crowding 
and saturation are not significant, and the completeness functions 
reflect this. The core of the cluster is imaged on the PC chip (solid 
line), and the completeness function for this chip matches well the 
completeness functions for the three WFC chips, which imaged the 
outer regions of the cluster. If incompleteness (or rapid spatial vari- 
ation of completeness) was a significant issue in the central regions, 
then the position averaged PC function would be significantly de- 
graded in comparison to the three WFC functions. The second clus- 
ter in Fig. 2, NGC 2005, is an example of this. The effect of this 
cluster's extremely compact core on the PC completeness function 
is clear, and it is evident that a positional resolution of 160 pix- 
els is not adequate to fiiUy describe the completeness variations. In 
fact, in the very central regions, the situation is even worse than 
this - the completeness is so low that the measurements become 
meaningless. To demonstrate this, from the three WFC functions, 
we expect the non-crowded regions of the PC to have c ~ 0.85 
at V555 ~ 24, whereas we measure a spatially averaged value of 
c ~ 0.45. From the PC image, approximately half of the chip may 
be non-crowded, implying that the integrated completeness in the 
crowded region must he, c ^ 0.05 at Vsss ~ 24 (cf. Fig. 3(b)). 
In a case such as this, we had to use additional information from 
a short exposure (see Section 4.2.3, below), otherwise the surface 
brightness profile became useless. 

We note that none of the completeness functions plotted in 
Fig. 2 ever reaches c = 1 - this is due to the integration over 
colour and position, since bad pixels and bright stars prevent all 
fake stars being recovered. Certainly however, individual bins in 
position-brightness-colour space can (and do) have full complete- 
ness. Completeness functions can also tell us a small amount about 
the background field near a cluster. The degradation evident in the 
WFC completeness functions for NGC 2005 is primarily due to the 
dense field population in the region of this cluster (NGC 2005 lies 
very near the LMC bar). By comparison, the field population near 
NGC 2213 is quite sparse. 

4.2.3 Short exposures 

As discussed in Section 4.2.2, in nine cases where severe crowding 
and/or saturation significantly degraded the completeness in parts 
of an image, we used information from short exposure images to 
complete the surface brightness profiles. This occurred specifically 
for the very old, compact clusters NGC 1754, NGC 1835, NGC 



1916, NGC 2005, and NGC 2019, and the very young clusters NGC 
1984, NGC 2011, NGC 2156, and R136. NGC 1786 required cor- 
rection but no additional data was available. 

To add the short exposure data is a relatively straightforward 
process. For a given cluster, the short exposure image is reduced 
in exactly the same manner as the original image, and a surface 
brightness profile constructed. When plotted over the original sur- 
face brightness profile, the short exposure profile is found to have 
a shape matching that of the original at large radii (where neither 
crowding nor saturation are an issue), but with a larger scatter and 
brighter by some tenths of a magnitude. This offset is a product 
of the shorter image having a brighter saturation limit and there- 
fore including more bright stars in its profile than the longer image, 
whereas the larger scatter is due to both the inclusion of many of 
the brightest stars in the cluster (see also Section 4.2.4) and the ex- 
clusion of large numbers of faint stars, which do not get measured 
on the short image. To eliminate the offset, we impose a brightness 
limit Viim on the short exposure photometry to cut out the brightest 
stars, and vary Viim until the profiles overlap in their outer regions. 
This is equivalent to removing the brightest stars in the cluster from 
the calculations, as discussed below in Section 4.2.4, and so also 
reduces the noise in the short exposure profile. It is then simple to 
measure the radius rdcv at which incompleteness becomes an issue 
and the two profiles begin to deviate, and form a composite profile 
using photometry for stars with r < rdev from the short exposure 
(with Viirn imposed) and the remainder (with r > rdev) from the 
long exposure. This keeps scatter to a minimum while alleviating 
the incompleteness at small radii. 

Fig. 3 shows an example of this process for NGC 2005. In 
the previous Section (4.2.2) we showed that the long exposure 
photometry for this cluster suffered serious incompleteness in the 
central region. Further evidence of this can be seen in Fig. 3(a). 
The two surface brightness profiles show a similar shape up to 
log rdeu ~ 0.8, or rdev ~ 6'.'3. At this point incompleteness in 
the long exposure photometry becomes significant and the surface 
brightness profile eventually turns over because of the amount of 
missing flux. Not only is the completeness very low in this region, 
but it must vary on a scale shorter than ~ 6", or ~ 130 pixels. To 
accotmt for the crowding, measurements are taken from the short 
photometry within 6'.'3. The advantage gained from this is evident 
in Fig. 3(b), which shows the completeness functions for the long 
and short exposure photometry, for r < rdev ■ As before, the com- 
pleteness functions have been integrated over position and colour 
to be a function of brightness only. The terrible incompleteness in 
the long exposure photometry is very evident. It is also clear that 
the short exposure photometry within the central region accounts 
very well for this incompleteness, and we are therefore justified in 
using only the short exposure photometry for the surface brightness 
profile in this region. 

4. 2. 4 Saturated stars 

At this stage we comment briefly on the treatment of saturated stars 
in our images. Such stars were not measured by multiphot and were 
therefore not included in the construction of surface brightness pro- 
files. For most clusters, saturated stars are not a significant factor. 
Of the 53 clusters in the sample, thirty had less than five saturated 
stars within approximately two core radii of the centre. When data 
from short exposures were included to circumvent extreme crowd- 
ing and saturation (see Section 4.2.3), this number rose to 39. For 
all of these clusters, the presence of so few saturated stars meant 
that their neglect did not significantly alter the surface brightness 
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Table 3. Clusters with more than five saturated stars within two core radii. 



Cluster 


logr 








(yr) 






NGC1466 


10.10 


15 


0.002 


NGC171 1 


7.70 


45 


0.018 


NGC1786 


10.18 


~ 20^^ 


■~ 0.008 


NGC18()5 


7.00 


30 


0.034 


NGC1818 


7.40 


35 


0.016 


NGC1841 


10.09 


~ 50'' 


0.003 


NGC1866 


8.12 


30 


0.009 


NGC1898 


10.15 


35 


0.006 


NGC2031 


8.20 


35 


0.012 


NGC2159 


7.60 


25 


0.022 


NGC2164 


7.70 


30 


0.016 


NGC2172 


7.60 


30 


0.033 


NGC2210 


10.20 


30 


0.005 


NGC2257 


10.20 


10 


0.001 



" Approximate number of saturated stars within two core radii. 

Ns divided by the total number of stars within two core radii. 
° This figure is an estimate since NGC 1786 suffers from extreme 
crowding and no short exposure was available. 

NGC 1841 has an extremely large core radius. 



profiles derived for these clusters, especially given that on average 
several thousand stars were measured in the central regions of each. 

The remaining fourteen clusters are listed in Table 3, along 
with an estimate of the total number of saturated stars within two 
core radii for each cluster. We felt comfortable leaving these stars 
out of the calculations, for the following reasons. Firstly, their pres- 
ence did not degrade the quality of the images or the measurements 
made from these images (otherwise short exposure data would have 
been obtained - see Section 4.2.3). Secondly, it is clear that for each 
cluster, the fraction listed in Table 3 is tiny (although larger than 
for the other 39 clusters). The fourteen clusters in Table 3 are either 
young (t ^ 100 M5T) or very old (t > 10 Gyr). In such clus- 
ters, the brightest stars are either giants or upper main-sequence 
stars and are in very brief (but very luminous) phases of evolution, 
especially compared to the dynamical timescales (e.g., the median 
relaxation time) of the clusters they are in. They therefore may not 
in general represent the spatial distribution of the underlying stellar 
population. It does not make sense to measure a surface brightness 
profile which is dominated by the output of an essentially randomly 
distributed tiny fraction of the stars in a cluster - in practice, this 
simply adds large amounts of noise to the profile. 

Removing these stars from the calculations provides less noisy 
profiles without compromising the measurement of structural pa- 
rameters. Elson and collaborators for example (Elson et al. (1989); 
E91; Elson (1992)) routinely remove the brightest stars when cal- 
culating their surface brightness profiles, in order to make sure the 
profile reflects the underlying stellar distribution. Specifically, El- 
son et al. (1989) found that the removal of these stars (in the con- 
struction of their "cleaned" profiles) reduced the scatter in their pa- 
rameter measurements by ^ 40 per cent on average without alter- 
ing the parameters by more than 5 per cent. Elson (1992) removed 
the brightest 120 stars per cluster in the construction of her cleaned 
profiles - again, this provided much smoother profiles without al- 
tering the derived structural parameters by a significant amount. 
Some experimentation on our part showed similar results. While 
we do not find it necessary to remove quite so many stars, in the 



cases listed in Table 3, it is evident from the final results presented 
in Fig. 6 that leaving the brightest (saturated) stars out of the con- 
struction process has resulted in clean, low noise profiles, which 
should reflect the underlying stellar distribution of each cluster. 

4.2.5 Error estimation 

To estimate the internal error ai for an annulus we divided it into 
eight sectors of equal area, and calculated the surface brightness in 
each using Eq. 1 with the factors Ai set to unity. The internal error 
for the aimulus is the standard deviation of the surface brightness 
values for these eight sectors (Djorgovski 1987). This technique en- 
sures that the internal error of an annulus reflects the error due to 
the discreteness of the light distribution - that is, clumps of stars or 
single bright stars will increase the surface brightness of an annu- 
lus, but also its error by an appropriate amount. This is essential in 
the case of LMC clusters, many of which are patchy and irregular. 
Similarly, the error for an annulus reflects the width of that annulus 

- important given our four different annulus sets. 

In the outer regions of most clusters however, we foimd that 
the error bars so calculated were considerably smaller than the 
random scatter in the points. As an example see Fig. 5, which 
demonstrates the background subtraction technique described in 
Section 4.3.2 but also well illustrates the present issue. The under- 
estimation is caused because Ai 3 for the outer annuli, so that 
while the flux for such an annulus is scaled to that for the full area, 
the errors are calculated over only a third of the full area and there- 
fore do not reflect any large scale variations for that annulus. This 
effect is coupled with field star contamination for many clusters. 
This background is negligible in the inner regions of a cluster, but it 
can become significant in the outer regions. The background tends 
to be much more smoothly distributed than the cluster itself, and so 
artificially shrinks the error bars determined using sectors. To solve 
this problem, for every annulus we also calculated the error from 
Poisson statistics. When this became significantly larger than the 
sector error, we took the Poisson error instead. This condition only 
ever became true in the outer regions of clusters, meaning that in 
the central regions of a cluster the errors are sector errors, while in 
the outer regions they are Poisson errors, with a smooth transition 
across some intermediate region (usually logr ~ 1.6). The Pois- 
son errors tend to over-estimate the scatter between points - this is 
due to the relatively low numbers of stars in the large outer annuli 

- but we find this preferable to an imder-estimate, especially given 
our process of background subtraction (see Section 4.3.2). 

Finally, we note that for some clusters, the (inner) sector errors 
are also larger than the RMS point-to-point scatter. In considering 
this, it is important that only the points in a single atmulus set are 
compared - consecutive points from different annulus sets use at 
least some of the same information, so one would expect their scat- 
ter to be small. One should also examine the brightness range of a 
profile - some cover up to 9 magnitudes while others only 2 or 3 
magnitudes. This can give the impression of enormous error bars 
on some profiles and tiny ones on others. Nonetheless, particularly 
in low surface brightness clusters such as NGC 1841, NGC 2121, 
NGC 2209, SL 663, and SL855, the effect persists. It is due to the 
low density of these clusters - the stars (especially the bright stars) 
are sparsely distributed. Therefore, while the surface brightness for 
each annulus is the azimuthal average over a comparatively large 
area, the eight sectors used to calculate the error for an annulus 
may have a considerable standard deviation in their surface bright- 
nesses and correspondingly result in errors which are larger than 
the RMS scatter between points. In the richer clusters (e.g., NGC 
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1466, NGC 1805, NGC 1831, NGC 2210, etc) the errors match the 
RMS point-to-point scatter well. 



4.3 Profile fitting and background subtraction 

4.3. 1 King versus EFF profiles 

Traditionally, the surface brightness profiles of old globular clusters 
are described by families of King-type models - for example the 
empirical King (1962) models: 

/^(r) = k \ i I (2) 

where n is the tidal radius of the cluster, and rc the core ra- 
dius. Provided rt » Vc, the core radius may be taken as the 
radius at which the surface brightness has dropped to half its cen- 
tral value. It is therefore useful to define the concentration param- 
eter c = log(rt/rc). For LMC clusters, concentration parame- 
ters of 1.0 ^ c ^ 2.0 are measured (cf. EFF87) implying ratios 
10 <, n/rc J$ 100. 

However, studies of young LMC clusters by Elson and col- 
laborators (EFF87; Elson et al. (1989); E91) have shown that these 
clusters do not appear to be tidally truncated, even at radii of several 
hundred arcseconds. EFF87 argue that the lack of tidal truncation in 
a young cluster is due to the cluster expanding over its Roche limit 
as a result of mass loss or violent relaxation. Given that a young 
cluster has only completed a fraction of its orbit about the LMC, 
any stars outside the Roche limit will not have been tidally stripped 
and instead surround the cluster in an unbound halo. In this case, a 
more suitable profile is: 

M(r)=Mo(^l + ^^ ' (3) 

where no is the central surface brightness, a is a measure of the core 
radius and 7 is the power-law slope at large radii. For 7 ~ 2 this is 
essentially Eq. 2 with rt 00. Wc hereafter refer to profiles of the 
form of Eq. 3 as EFF models, after Elson, Fall & Freeman (1987), 
who first introduced them into star cluster research. Fig. 4 shows an 
empirical King model compared with an EFF model, both having 
been calculated with a core radius of rc = 10", typical of the LMC 
clusters measured by EFF87 and Mateo (1987). For the King profile 
we take c = 1.5. It is evident that there is little deviation between 
the two profiles in the iimer regions. The difference only becomes 
large at radii approaching rt. 

The presence of both very young and very old clusters in our 
sample raises a potential dilemma over the fitting of one type of 
profile or the other, particularly for intermediate age clusters. For- 
tunately this is solved for us by the small field of view of WFPC2, 
which limits our surface brightness profiles to have maximum radii 
of r ~ 80" or log r ~ 1.9. From Fig. 4 we see that at such radii 
the two profiles are essentially indistinguishable, and so we may fit 
either. Because of the virtual impossibility of obtaining any accu- 
rate information about rt without measuring to at least this radius, 
we chose to fit EFF profiles to our data. 

The parameter a in an EFF profile may be linked to the King 
core radius by setting fj, = /lio/2 in Eq. 3 and rearranging to obtain 

re = ai2'/^ - 1)1/^ (4) 

We note that this assumes n >> rc, which is valid for most of the 
clusters in our sample. For several however (e.g., NGC 1841 and 



NGC 2257) rc is very large and 7 quite steep. In such cases, Eq. 4 
may imderestimate the core radius by 3 to 10 per cent (Elson 1992). 

4.3.2 Fitting procedure and background subtraction 

When fitting the profiles described in the previous section, in most 
cases it was necessary to correct for field star contamination. Espe- 
cially for clusters near the LMC bar, the stellar background contri- 
bution is not negligible and tends to flatten out the cluster profiles 
in their outer regions (see e.g.. Fig. 3). Because 7 is largely de- 
termined from the slope of a profile well away from the core and 
is therefore very sensitive to field star contamination, we had to 
subtract the background level in order to avoid large systematic er- 
rors. Traditionally one would use an observation well away from 
the cluster to obtain the local background, and subtract this level 
from each annulus. However, we had no access to such background 
fields for our entire sample, nor were any of our fields large enough 
in area to allow background determinations away from a cluster. 
We therefore had to try and determine the background levels indi- 
vidually from each surface brightness profile. 

One way of doing this is to fit Eq. 3 but with an extra parame- 
ter (j){r) added - we set as a constant to represent a simple back- 
ground. However, these four-parameter fits did not work as well 
as desired - they tended to over-estimate the background, causing 
7 to be too steep by factors of up to three. The reason for this is 
not clear; it would seem that the creation of an additional degree of 
freedom degraded the quality of the fit for the original three param- 
eters. We therefore implemented a recursive method using the fact 
that for r >> a, Eq. 3 with a background contribution becomes: 

« Mo (0 + (t> 

« ^o(i_)""(22/7_i)-7/2 + ^ (5) 

with Eq. 4 substituted for a. 

We first fit an EFF model to the inner region of a cluster pro- 
file using the two narrow annulus sets, to estimate no and rc. In 
this region (r ^ 30") the background contribution is negligible, 
in all cases being at least 5 magnitudes below the central surface 
brightness. The background therefore does not affect a measure- 
ment of the central surface brightness /io nor the core radius rc, 
which is essentially the shape of the profile in the inner region. 
Hence it is a good approximation to determine these without sub- 
tracting any background. We next use these estimated values and 
fit a two-parameter model (7, <j)) of the form of Eq. 5 to the outer 
region of the profile, using the two wider annulus sets and ensur- 
ing r » a. This allowed us to determine </> and subtract it from 
every annulus. After the background subtraction, we calculate the 
Poisson errors for each point and substitute as appropriate (see Sec- 
tion 4.2.5). Finally, we fit another EFF model to the new subtracted 
profiles, using all four annulus sets, and thereby determine the best- 
fitting parameters. This method allows us to treat severely contam- 
inated clusters (e.g., NGC 1898 and NGC 1916) equally with clus- 
ters that have no evident background contribution (e.g., Hodge 1 1 
and NGC 1841). 

To test the approximation described above, we examined the 
final parameters (no, rc) to see how well they corresponded to the 
determinations from the central fit. In no cases did we observe sig- 
nificant deviation, suggesting that the approximation is good. In 
addition, we checked for any systematic errors introduced by this 
method by comparing our determinations of 7, which is very sensi- 
tive to the background level, with previously pubUshed work. This 



© 2001 RAS, MNRAS 000, 000-^00 



Structure of rich LMC clusters 1 1 



is presented in Section 5.1.1, and suggests that we introduced no 
significant systematic errors. A final point concerns the effect of 
non-constant backgrounds such as might be encountered near the 
bar - in this case, cj) = Such a background would have two 
effects given the procedure described above. The first would be to 
cause a very small inaccuracy in our determined centre. The second 
would be to increase the noise in the outer region of a profile, be- 
cause of the fact that we image only a fraction of any outer annulus. 
This in turn would cause a larger random error in the determination 
of 7, but this would be accounted for by our bootstrapping method 
of determining errors in the derived parameters (see below). 

Fig. 5 shows an example of the background fitting, again for 
NGC 2005. In this case, the estimated parameters from the initial 
fit were a ~ 4" and fio ~ 16.9, and the background measured as 
Vbg = 21.83. As mentioned in Section 4.2.2, NGC 2005 is very 
near the LMC bar, and has a densely populated background field. 
The derived background intensities for this cluster and the others 
in the sample are consistent with the LMC isophote maps of de 
Vaucouleurs (1957). Equivalently, a plot of background intensity 
vs. Ropt shows the expected decrease with increasing R. 

In implementing each fit, we develop a grid in parameter space 
and choose the combination of parameters which minimizes the 
weighted sum: 

i=i \ / 

where /li is the surface brightness of the i-th annulus, ai is the 
error in this value, A'^^ is the total number of annuli in the set in 
question, and the other parameters are as defined in Eq. 3. A re- 
fined mesh is then expanded about this parameter combination and 
the iteration continued until convergence. Given that our parameter 
space is topologically well behaved, this method was accurate and 
efficient, with convergence typically within ten iterations. Because 
our technique of using four annulus widths effectively counts each 
star four times, to maintain independence between data points we 
fit each annulus set individually as appropriate. To obtain the best 
fit curve for an iteration we average the individual fits, weighted 
by their values. We determined errors in the three parameters 
obtained from the final fits to the background-subtracted profiles 
using a bootstrap method (Press et al. 1992, p691) with 1000 recur- 
sions per fit. Again, this was done separately for each annulus set to 
maintain the independence of the points. The errors so determined 
represent random errors in the best-fit parameters. Systematic er- 
rors such as those potentially introduced by the backgrotmd fitting 
are not included in the estimation. 



5 RESULTS 

5.1 Profiles and structural parameters 

The background- subtracted F555W surface brightness profiles for 
each of the 53 clusters are presented in Fig. 6. We plot each of 
the four annulus sets on the same axes, to demonstrate the high de- 
gree of consistency between them. For each cluster, the best fit EFF 
profile is also plotted, the core radius indicated, and the best fit pa- 
rameters listed. These results are sunmiarized in Table 4 along with 
their corresponding errors, the calculated centre of each cluster and 
the maximum radial extent Vm of each profile. Three of the pro- 
files (NGC 1754, NGC 1786, NGC 1916) are incomplete in their 
inner regions - this is due to the effects of crowding, even includ- 
ing supplementary photometry from short exposures (a short expo- 



sure was unavailable for NGC 1786, however). The derived core 
radii for these clusters, and for the clusters NGC 1835, NGC 2005, 
NGC 2019 and R136 which also required short exposure photome- 
try, should be considered upper limits. 

This is one of the largest published studies of LMC cluster sur- 
face brightness profiles, and as far as we are aware, the only one to 
use ffiTdata. In addition, approximately one quarter of the clusters 
in the sample do not have previously published surface brightness 
profiles. There is much to be learned from such detailed data, cov- 
ering many aspects of globular cluster and LMC astronomy. We 
can however split the sample into sub-groups and obtain some im- 
mediate and interesting results, which we discuss below. We first 
compare our results with those from the several other published 
large-scale studies. 

5.1.1 Comparison with previous woric 

At least three-quarters of the clusters studied in this paper have 
published surface brightness profiles, concentrated mainly in four 
large studies - those of EFF87, Mateo (1987), E91 (see also El- 
son et al. (1989)), and Kontizas and collaborators (Kontizas et al. 
1987a; Kontizas et al. 1987b; Metaxa et al. 1988; Chrysovergis et 
al. 1989). It is useful to compare our results to those from these pa- 
pers and thereby establish the validity of our reduction procedure. 

In particular, our background subtraction procedure might in- 
troduce systematic errors into the measured structural parameters. 
We checked this by comparing our derived values for 7 with those 
from EFF87 for the 10 young clusters (r < 3 x 10* yr) in their 
sample, since 7 is very sensitive to the subtracted background. The 
comparison is plotted in Fig. 7. The profiles from EFF87 extend 
to typically r,„ ^ 250" and their background measurements are 
made at these large radii. The agreement with our fitted values of 7 
is generally good, and within the errors for most points. If anything 
there is a slight tendency for our 7 measurements to be slightly 
larger than those from EFF87 (i.e., we slightly over-estimated the 
background - as would be expected), but this discrepancy is not 
significant. EFF87 find a range of 2.2 <, ^ <, 3.2 and a me- 
dian value 7 w 2.6. For the same ten clusters we find a range of 
2.26 7 3.45 but a median value 7 « 2.95. We note however 
that the sample size is barely large enough for meaningful statis- 
tics. Extending the average to cover the 22 clusters in our sample 
with T < 3 X 10** yr we obtain a range 2.01 7 ^ 3.79 and a 
median 7 w 2.59. We are therefore confident that our backgrotmd 
subtraction method has not introduced any significant systematic 
errors into the structural parameter measurements. 

We have also taken the core radii from the first three stud- 
ies listed above for all the clusters in corrmion with our sample, 
and plotted these against our measured values (Fig. 8). Again, the 
agreement is relatively good; however there is some considerable 
scatter in the plot. Nonetheless, there are good explanations for 
many of the errant points. Our core radii are generally somewhat 
larger than those measured by EFF87 (filled triangles) - however 
the irmer regions of the profiles from EFF87 are a combination of 
literature data and photomultiplier aperture measurements and the 
errors quoted for a are up to 50 per cent and more, which would 
bring most if not all of our results well into the error margins. 
More encouraging is the good agreement of our measurements with 
those from E91 (filled circles), again for relatively young clusters 
(10^ < r < 10" yr). The profiles from this study were measured 
from CCD images, with backgrounds determined at r ~ 200". 
The most errant point (falling well below the equality line) is for 
the young cluster NGC 2100. We speculate that the different satu- 
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Table 4. Structural parameters for the cluster sample derived from the best-fitting F555W EFF profiles. 



Cluster 

Name 


Centre (J2000.0)" 

a S 


/U655 (0)' 


b 


a 

(") 






7 






r-c 

(") 




(pc)° 




(") 


NGC1466 


03'"44'"32!9 


-71°40'13" 


19, 


,08 ± 0, 


,05 


14.89 ± 0.79 


3, 


,31 


± 


0. 


12 


10.73 ±0.36 


2, 


,61 ± 0, 


,09 


76 


NGC1651 


04'i37™31!l 


-70°35'02" 


20, 


,38 ± 0, 


,04 


15.82 ± 1.26 


2, 


,21 


± 


0. 


12 


14.75 ±0.70 


3, 


.58 ± 0, 


,17 


76 


NGC1711 


04'' 50"" 37! 3 


-69°59'04" 


18, 


,28 ± 0, 


,05 


10.93 ± 0.83 


2, 


,78 


± 


0. 


13 


8.78 ± 0.44 


2, 


.13 ±0, 


,11 


72 


NGC1718 


04''52'"25?6 


-67°03'06" 


19, 


,18 ±0, 


,07 


9.38 ± 0.76 


2, 


.31 


± 


0. 


10 


8.52 ± 0.49 


2, 


.07 ±0. 


,12 


76 


NGC1754 


04''54'"18?9 


-70°26'31" 


17, 


,48 ± 0, 


,26 


4.12 ±0.65 


2, 


.43 


± 


0. 


,09 


3.61 ± 0.49 





.88 ± 0. 


,12 


76 


NGC1777 


04/1 55m, 4gsg 


-74° 17' 03" 


20, 


,01 ± 0, 


,02 


16.62 ± 0.89 


2, 


.60 


± 


0. 


,10 


13.96 ± 0.44 


3 


.39 ± 0, 


,11 


76 


NGC1786 


04''59'"07?9 


-67° 44' 45" 


17, 


,46 ± 0, 


,11 


6.79 ± 0.67 


2, 


,78 


± 


0. 


13 


5.45 ± 0.38 


1, 


.33 ± 0, 


,09 


76 


NGC1805 


05''02'"21!8 


-66°06'42" 


18, 


,01 ±0, 


,06 


6.84 ± 0.42 


2, 


,81 


± 


0, 


,10 


5.47 ±0.23 


1 


.33 ± 0, 


,06 


69 


NGC1818 


05^04'" 13!8 


-66°26'02" 


18, 


,36 ± 0, 


,05 


12.50 ± 0.78 


2, 


,76 


± 


0, 


12 


10.10 ±0.39 


2 


.45 ± 0, 


,09 


76 


NGC1831 


05''06'"17!4 


-64°55'11" 


19, 


,08 ± 0, 


,04 


25.81 ± 1.59 


3, 


,41 


± 


0, 


19 


18.28 ±0.59 


4, 


.44 ± 0, 


,14 


76 


NGC1835 


05''05'"06!7 


-69°24'15" 


16, 


,37 ± 0, 


,08 


6.35 ±0.46 


3, 


,11 


± 


0, 


10 


4.76 ±0.26 


1 


.16 ± 0, 


,06 


76 


NGC1841 


04''45'"23!9 


-83°59'56" 


21, 


,59 ± 0, 


,02 


53.62 ± 5.07 


4, 


,55 


± 


0. 


,61 


32.00 ± 0.72 


7, 


.77 ± 0, 


,17 


78 


NGC1847 


05''07'"07?7 


-68°58'17" 


18, 


,54 ± 0, 


,06 


8.58 ± 0.44 


2, 


.05 


± 


0. 


,13 


8.44 ± 0.34 


2, 


.05 ± 0, 


,08 


42 


NGC1850 


05''08'"41?2 


-68°45'31" 


16, 


,70 ± 0, 


,07 


11.11 ± 1.03 


2, 


.18 


± 


0. 


,11 


10.48 ± 0.64 


2 


.55 ±0, 


,16 


69 


NGC1856 


05''09"'31?5 


-69°07'46" 


16, 


,68 ± 0, 


,07 


6.81 ±0.46 


2, 


.01 


± 


0. 


,05 


6.80 ± 0.36 


1, 


.65 ± 0, 


,09 


76 


NGC1860 


05''10'"38!9 


-68°45'12" 


19, 


,72 ± 0, 


,11 


9.74 ± 0.53 


2, 


,25 


± 


0. 


13 


9.00 ± 0.50 


2, 


.19 ±0, 


,12 


76 


NGC1866 


05''13'"38!9 


-65°27'52" 


18, 


,52 ±0, 


,02 


18.60 ±0.72 


3, 


,04 


± 


0, 


,09 


14.15 ±0.32 


3, 


.44 ± 0, 


,08 


76 


NGC1868 


05'' 14™ 36! 2 


-63° 57' 14" 


17, 


,99 ± 0, 


,04 


8.83 ±0.41 


3, 


,07 


± 


0, 


10 


6.67 ±0.19 


1 


.62 ± 0, 


,05 


76 


NGC1898 


05''16'"42!4 


-69°39'25" 


19, 


,05 ± 0, 


,06 


8.79 ±0.91 


2, 


,14 


± 


0, 


14 


8.40 ±0.55 


2 


.04 ± 0, 


,13 


76 


NGC1916 


05''18'"37!5 


-69°24'25" 


15, 


,93 ± 0, 


,09 


4.08 ±0.28 


2, 


,70 


± 


0, 


,08 


3.35 ±0.18 





.81 ± 0, 


,04 


76 


NGC1984 


05^27'" 40!8 


-69°08'05" 


18, 


,62 ± 0, 


,11 


4.44 ± 0.41 


2, 


,27 


± 


0. 


11 


4.07 ±0.28 





.99 ± 0, 


,07 


76 


NGC2004 


05''30'"40!9 


-67° 17' 09" 


16, 


,93 ± 0, 


,08 


7.56 ±0.86 


2, 


,53 


± 


0. 


,12 


6.47 ± 0.53 


1, 


.57 ±0, 


,13 


76 


NGC2005 


05''30'"10!3 


-69°45'09" 


17, 


,04 ± 0, 


,13 


4.33 ± 0.54 


2, 


.60 


± 


0. 


,14 


3.63 ± 0.33 





.88 ±0, 


,08 


69 


NGC2011 


05''32'"19?6 


-67°31'14" 


18, 


,93 ± 0, 


,15 


5.17 ±0.60 


2, 


.22 


± 


0. 


09 


4.81 ± 0.51 


1, 


.17 ±0, 


,12 


76 


NGC2019 


05''31'"56?6 


-70° 09' 33" 


16, 


,73 ±0, 


,10 


4.21 ±0.37 


2, 


.52 


± 


0. 


,09 


3.61 ± 0.24 





.88 ± 0, 


,06 


76 


NGC2031 


05''33'"41!1 


-70°59'13" 


18, 


,93 ± 0, 


,06 


11.13 ±0.98 


2, 


,09 


± 


0. 


11 


10.81 ±0.61 


2, 


.63 ± 0, 


,15 


76 


NGC2100 


05''42'"08!6 


-69°12'44" 


16, 


,26 ±0, 


,14 


5.73 ± 0.75 


2, 


,44 


± 


0, 


14 


5.02 ±0.50 


1 


.22 ± 0, 


,12 


76 


NGC2121 


05''48™11!6 


-71°28'51" 


20, 


,90 ± 0, 


,02 


33.83 ±3.83 


2, 


,25 


± 


0, 


,28 


31.22 ± 1.14 


7 


.59 ±0, 


,28 


76 


NGC2136 


05''52'"58!0 


-69°29'36" 


17, 


,37 ±0, 


,09 


12.33 ± 1.06 


3, 


,79 


± 


0, 


,23 


8.19 ±0.43 


1 


.99 ± 0, 


,10 


69 


NGC2153 


05''57™51!2 


-66°23'58" 


19, 


,17 ±0, 


,09 


4.71 ±0.79 


2, 


,64 


± 


0, 


,22 


3.91 ±0.44 





.95 ± 0, 


,11 


69 


NGC2155 


05''58'"33!3 


-65°28'35" 


20, 


,44 ± 0, 


,05 


20.30 ±2.20 


2, 


,85 


± 


0, 


,25 


16.06 ±0.98 


3 


.90 ±0, 


,24 


76 


NGC2156 


05''57'"49!7 


-68°27'41" 


18, 


,98 ± 0, 


,07 


8.15 ±0.51 


2, 


,79 


± 


0. 


,09 


6.53 ± 0.31 


1, 


.59 ± 0, 


,08 


76 


NGC2157 


05''57'"32?4 


-69°11'49" 


17, 


,78 ±0, 


,05 


13.76 ±0.91 


3, 


.45 


± 


0. 


,18 


9.68 ± 0.36 


2 


.35 ± 0, 


,09 


76 


NGC2159 


05''58"'03?0 


-68°37'27" 


19, 


,66 ± 0, 


,05 


12.73 ±0.94 


3, 


.00 


± 


0. 


,15 


9.77 ± 0.47 


2 


.37 ± 0, 


,11 


69 


NGC2162 


06''00'"30?4 


-63°43'19" 


19, 


,92 ± 0, 


,08 


12.97 ± 1.40 


2, 


,91 


± 


0. 


23 


10.13 ±0.66 


2 


.46 ± 0, 


,16 


76 


NGC2164 


05''58'"55!9 


-68°31'00" 


18, 


,22 ± 0, 


,05 


10.16 ±0.43 


2, 


,96 


± 


0, 


,07 


7.85 ±0.24 


1 


.91 ± 0, 


,06 


72 


NGC2172 


OghoO™06!4 


-68°38'15" 


20, 


,21 ± 0, 


,07 


13.77 ± 0.97 


2, 


,94 


± 


0, 


14 


10.68 ±0.49 


2 


.60 ± 0, 


,12 


76 


NGC2173 


05''57'"58!5 


-72°58'40" 


19, 


,86 ±0, 


,05 


13.13 ± 1.35 


2, 


,37 


± 


0, 


16 


11.70 ±0.74 


2 


.84 ±0, 


,18 


76 


NGC2193 


06''06'"17!0 


-65°05'54" 


19, 


,77 ±0, 


,08 


7.65 ±0.71 


2, 


,50 


± 


0, 


11 


6.59 ±0.43 


1 


.60 ± 0, 


,10 


72 


NGC2209 


06''08'"34!8 


-73° 50' 12" 


22, 


,09 ± 0, 


,06 


22.97 ±2.39 


2, 


,08 


± 


0, 


16 


22.34 ± 1.34 


5 


.43 ± 0, 


,33 


76 


NGC2210 


06''11'"31!5 


-69°07'17" 


18, 


,18 ±0, 


,05 


11.75 ±0.56 


3, 


,51 


± 


0. 


,12 


8.18 ± 0.24 


1, 


.99 ± 0, 


,06 


76 


NGC2213 


nfi'»i n'"42^2 

\J\J Xl_/ '-^Zi* 


— 71°31'46" 


19, 


,16 ± 0, 


,05 


9.36 ± 0.69 


2, 


.84 


± 


0. 


,13 


7.43 ± 0.36 


1, 


.80 ± 0, 


,09 


76 


NGC2214 


06''12'"55?8 


-68° 15' 38" 


18, 


.36 ± 0, 


,07 


9.55 ± 0.99 


2, 


.26 


± 


0. 


,14 


8.79 ± 0.59 


2 


.14 ±0, 


,14 


76 


NGC2231 


O6''20'"42?7 


-67°31'10" 


20, 


,48 ± 0. 


,12 


11.42 ±2.79 


2, 


,07 


± 


0. 


26 


11.14 ± 1.55 


2 


.71 ±0, 


,38 


76 


NGC2249 


06''25'"49!8 


-68°55'13" 


19, 


.05 ± 0, 


,04 


12.13 ± 0.83 


3 


,40 




0, 


17 


8.61 ± 0.35 


2 


.09 ± 0, 


,09 


76 


NGC2257 


06 ''30™ 12! 1 


-64°19'42" 


20, 


,93 ± 0, 


,03 


38.61 ± 2.71 


3, 


,54 


± 


0, 


,27 


26.75 ±0.73 


6 


.50 ±0, 


,18 


76 


SL663 


05''42™28!8 


-65°21'44" 


22, 


,40 ± 0, 


,05 


28.83 ±2.24 


2, 


,13 


± 


0, 


15 


27.60 ± 1.27 


6 


.71 ±0, 


,31 


76 


SL842 


06''08"'14!9 


-62° 59' 15" 


21, 


,07 ± 0, 


,11 


14.55 ± 1.75 


3, 


,37 


± 


0, 


,29 


10.37 ±0.76 


2 


.52 ± 0, 


,18 


69 


SL855 


06'' 10™ 53! 7 


-65°02'30" 


22, 


,93 ± 0, 


,03 


23.10 ±2.91 


2, 


,12 


± 


0, 


,23 


22.18 ± 1.35 


5 


.39 ±0, 


,33 


76 


H0DGE4 


05''32"'25!2 


-64°44'11" 


20, 


,59 ± 0, 


,06 


15.84 ±1.65 


2, 


.13 


± 


0. 


,15 


15.20 ± 1.08 


3 


.69 ± 0, 


,26 


76 


HODGEll 


06''14"'22?3 


-69°50'50" 


19, 


,60 ± 0, 


,05 


13.66 ±1.17 


2, 


.38 


± 


0. 


,12 


12.14 ±0.66 


2 


.95 ± 0, 


,16 


76 


H0DGE14 


05/t 28m 3933 


-73°37'49" 


20, 


,28 ±0, 


,09 


8.38 ±0.93 


2, 


.41 


± 


0. 


,16 


7.39 ± 0.56 


1, 


.80 ± 0, 


,14 


76 


R136 


05''38'"42!5 


-69°06'03" 


12, 


,66 ± 0, 


,09 


1.48 ±0.12 


2, 


,43 


± 


0. 


09 


1.30 ±0.08 





.32 ± 0, 


,02 


10 



" We find our centering algorithm to be repeatable to approximately ±1", notwithstanding image header inaccuracies (see Section 4.1). Given this 
precision, coordinates in 5 are provided to the nearest arcsecond. Those in a are reported to the nearest tenth of a second, but the reader should 
bear in mind that at 5 = —69°, one second of RA corresponds to approximately five seconds of arc - in other words, the uncertainty in a is 
approximately ±0?2. 
The V555 magnitude of one square arcsecond at the centre of a given cluster. 

When converting to parsecs we assume an LMC distance modulus of 18.5 which equates to a scale of 4.116 arcsec pc~^. 
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ration limits between the two studies (12^V^19 from E91 as 
compared witii 15 ^ V555 ^ 22 from the present study) might be 
the root of this difference. Furthermore, we have adopted the core 
radius measurements from E91 not corrected for seeing, because 
there is some evidence of over-correction in the quoted values. The 
seeing corrections move the errant point considerably to the left on 
the plot, but also move most of the other points to the left and away 
from the equaUty line. 

Both the studies by EFF87 and E91 describe the profiles of 
many young LMC clusters as being very irregular in comparison to 
the smooth profiles typical of old globular clusters. We also observe 
such irregularities, in the same forms as those noted by E91 - that 
is, bumps (e.g., NGC 2004), steps and wiggles (e.g., NGC 1711, 
NGC 1856), sharp shoulders or dips near the core radius (e.g., NGC 
1850, NGC 1860, NGC 1984) and central dips (e.g., NGC 2031). 

The measurements from the study of Mateo (1987) for old 
LMC clusters (filled squares), are from S-band images and so are 
not strictly comparable with ours; however we proceed for the sake 
of completeness. Again, there are no significant systematic differ- 
ences, and most errant points can be explained. The four smallest 
core radii from Mateo all correspond to upper limit measurements 
in our study - so the points likely lie closer to the equality line than 
indicated. The two errant points to the right of the plot are estimates 
taken from the literature - the exact method of estimation is not 
clear, but we conjecture that it might be a significant cause of the 
observed differences. We also note good agreement with the core 
radius for NGC 2257 (25'.' 1 as compared with our measurement of 
26'.' 75 ± 0'.'73), which falls outside the plot. 

It is more difficult to reconcile our results with those from 
the four studies of Kontizas et al. (1987a; 1987b), Metaxa et al. 
(1988), and Chrysovergis et al. (1989), which use density profiles 
from number counts. While this again means the derived parame- 
ters arc not strictly comparable with ours, the two sets should be 
roughly similar. However, we observe a large, apparently random 
scatter between them. The authors report that their core radius es- 
timates are derived from King model fits to the density profiles, 
which yield rt and the concentration c, and that the errors are ~ 15 
per cent in rt and ±0.25 in c, which is ~ 20 per cent in the typical 
value c w 1.5. We estimate that these errors could cause uncer- 
tainties of the order of ±60 per cent or more in derived values of 
rc, consistent with the difficulties associated with counting stars in 
crowded regions on photographic plates and the resolution of the 
profiles (most do not extend within ~ 25"). Such uncertainties can 
accoimt for most of the observed scatter. Nonetheless, several mea- 
surements are discrepant by factors of up to 10 (e.g., rc « 55" for 
NGC 1805 as compared with our measurement of rc ~ 5'.' 5) and 
we are unable offer a plausible explanation to fully reconcile these. 

5. 1.2 Double clusters and profile bumps 

There is a small amount of evidence for double clusters in our sam- 
ple. In particular, NGC 1850 is a well known double cluster, with 
a much smaller and younger companion - NGC 1850B - separated 
by approximately 32" (see e.g., Gilmozzi et al. (1994)). It is not 
clear that NGC 1850B is physically associated with NGC 1850 - it 
may be a chance superposition of two clusters along the same line 
of sight. Measurements of the mass of each cluster, and their three 
dimensional separation and relative velocity are required to settle 
this issue. 

Nonetheless, NGC 1850B appears in our images, and its pres- 
ence can be seen in the surface brightness profile for NGC 1850 as 
a bump at log r ~ 1.55. We also show our WFC2 l/-band image of 



NGC 1 850 in Fig. 9 - in this image, the centre of NGC 1 850 is to the 
lower left, and NGC1850B is evident to the upper right of the main 
cluster. The over-sampUng of our data has the effect of smearing 
out sharp features in the surface brightness profiles. For example, 
a very bright star will appear once in each annulus sot, resulting 
in four points at slightly different radii and giving the impression 
of a bump rather than a spike. This makes the bimip in the profile 
of NGC 1850 almost indistinguishable from that which would be 
caused by a bright star, although upon examination of the actual 
image, the nature of NGC 1850B is immediately apparent. Our un- 
published CMD for NGC 1850 also shows evidence of this young 
sub-cluster in the form of a few main sequence stars above the turn- 
off for the main cluster. We also note that because each annulus is 
essentially an azimuthal average, the amplitude of the bump due 
to NGC 1850B is reduced - that is, the peak surface brightness of 
NGC 1850B is certainly much greater than indicated in the surface 
brightness profile for NGC 1850. 

The above result prompted us to search the other cluster im- 
ages and profiles (correlating the two against each other) for similar 
evidence of double clusters, even though only two of our sample, 
NGC 2011 and NGC 2136, appear in the catalogue of LMC dou- 
ble clusters published by Bhatia et al. (1991). None of the images 
showed any evident double clusters. Surprisingly, we did however 
locate two very prominent bumps in the profiles of NGC 2153 and 
NGC 2213, at logr ~ 1.15 (14") and logr ~ 1.25 (18") re- 
spectively. We show DSS2 ii-band images of each of these clusters 
in Fig. 9. It is clear from these images, and from our WFPC2 im- 
ages, that secondary clusters are not the cause of these bumps. Fur- 
thermore, although it is evident from the images that each cluster 
contains bright stars, neither of the bumps can be due to a single 
bright star, because each bump appears in at least two consecutive 
points from single annulus sets (the l'.'5, 2" and 3" sets for NGC 
2156, and the l'.'5 and 2" sets for NGC 2213). It is therefore not 
clear what the cause of either bump is. It is however clear that the 
binarity of a cluster cannot be judged from its surface brightness 
profile, because these examples show that even single clusters can 
have significant bumps in their profiles. 

One possibility which is evident from the DSS images (but not 
necessarily so clear on the high resolution WFPC2 images) is that 
the chance positioning of several bright stars at exactly the right 
distance from a given cluster's centre but with different position 
angles could result in a bump. For NGC 2153, there is an arc of 
stars to the north of the main cluster and a fainter arc to the south 
which are at approximately the right distance to match the bimip 
in the profile. There is also a bright star to the south, but this is 
saturated on our WFPC2 images. Whether these arcs are simply 
random placement of several stars, or a physical association is not 
clear. For NGC 2213, it would appear that there are several bright 
stars at approximately the right distance from the centre but at dif- 
ferent angles about the cluster, which might cause the bump in its 
profile. Again, whether any physical significance can be attached 
to this bump is uncertain. 

In addition, three other clusters - NGC 2004, NGC 2155, and 
NGC 2159 - have less statistically significant bumps in their pro- 
files, at logr ~ 1.25 (18"), logr - 1.3 (20"), and logr ~ 1.3 
(20") respectively. We show DSS2 i?-band images of these clus- 
ters, also in Fig. 9. Just as for NGC 2153 and NGC 2213, it is not 
immediately clear what the causes of these bumps are. Certainly 
they are not due to single bright stars, and again there are no clear 
sub-clusters like NGC 1850B. NGC 2004 shows a ring of bright 
stars about its centre, but these are at too great a radius to cause the 
observed bimip. We do note a sUght extension of the main cluster 
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towards the upper left - this is also apparent in the WFPC2 image 
and might be the source of the bump in the profile. Again, the phys- 
ical significance of this extension is not clean For NGC 2155, there 
is no evident source for the bump, although this cluster does show 
a clumpy and dispersed structure. Finally, for NGC 2159, again the 
bump may be caused by several bright stars at the right distance 
from the centre but at different position angles. There is no other 
structure evident which might cause the bimip in its profile. 

The profiles for NGC 2011 and NGC 2136, the only clusters 
in our sample apart from NGC 1850 to appear in the catalogue by 
Bhatia et al. (1991) do not show any structure that might be due to 
their binary nature. However, both companions lie outside the re- 
spective fields of view for these clusters. NGC 1818 also likely has 
a small companion (see e.g., de Grijs et al. (2002a)) but the separa- 
tion is ~ 90" and it again lies outside our field of view. NGC 2214 
is regarded by some authors as a binary cluster in the late stages 
of a merger (Bhatia & MacGillivray 1988), and has been noted to 
have an anomalous core (Meylan & Djorgovski 1987). We do not 
observe any significant bumps in our surface brightness profile for 
this cluster, or any interesting structure in our image, but we do see 
a similar effect to that described by Meylan & Djorgovski where 
the surface brightness profile in the core appears similar to that for 
a post core-collapse cluster. That NGC 2214 is such a cluster is 
unlikely given its age (~ 60 Myr). Finally, NGC 2156 is rendered 
interesting because of the results of de Oliveira et al. (2000) who 
measure this cluster to have somewhat elliptical isophotes. Upon 
comparison with the results of iV-body simulations of star cluster 
encounters, they conclude that this ellipticity could result from a 
binary cluster interaction, given a favourable viewing angle. More- 
over, the young age of NGC 2156 combined with the behaviour 
of the simulation suggests that the merging cluster will not yet 
have been disrupted. However, we observe no evidence for such 
a distinct sub-cluster in our images, or any structure in the surface 
brightness profile. 

5.1.3 Post core-collapse clusters 

We find strong evidence for several post core-collapse (PCC) clus- 
ters amongst the old clusters in our sample. Theoretical studies, as 
well as observations of the Galactic globular cluster system have 
shown that the surface brightness profiles of PCC clusters are char- 
acterized by power-law cusps in their central regions, rather than 
ordinary King-type models. We have identified two clusters with 
profiles which clearly show such cusps as well as breaks to the 
power-law regions, and another two clusters which show cusps and 
breaks at lower signficance. A further three clusters have incom- 
plete profiles in their central regions (a result of severe crowding 
in these regions - see Section 4.2.2 and the beginning of Section 
5) but do appear to show breaks similar to the other four PCC can- 
didates. The three groups of clusters are, respectively: NGC 2005 
and NGC 2019; NGC 1835 and NGC 1898; and NGC 1754, NGC 
1786, and NGC 1916. 

We first address NGC 2005 and NGC 2019, which have pro- 
files clearly showing PCC-like cusps. This has previously been 
noted by Mateo (1987) for both clusters, and by Meylan & Djor- 
govski (1987) for NGC 2019. Mateo estimates power-law slopes of 
/3 ~ 1 for both clusters, where we assume the profile goes as 
in the central regions. He also notes sharp breaks in both profiles at 
the transition to the power-law region, which occur at 4'.'4±0'.'4 for 
NGC 2005 and 3'.'6 ± 0'.'4 for NGC 2019. We observe well defined 
power-law regions for both clusters, with sharp breaks at 5'.'6 ± 0'.'3 
for NGC 2005 and S'.'O ± 0'.'6 for NGC 2019, both slightly larger 



than the measurements of Mateo. If we fit power-law profiles to 
these inner sections, we obtain slopes of (5 = 0.74 for NGC 2005 
and /3 = 0.71 for NGC 2019. The profiles and power-law models 
are shown in Fig. 10(a) and (b), together with the best fitting EFF 
models. Clearly, for each cluster the power-law is a better fit in the 
centre than the EFF model. Further, it is interesting that studies of 
Galactic PCC clusters have measured slopes in the range 0.6 to 0.8 
(Lugger, Cohn & Grindlay 1995), or with a median value (3 w 0.9 
(Djorgovski & King 1986), and many of the PCC Galactic clusters 
also show breaks in their profiles (see Fig. 2 in Djorgovski & King, 
and Fig. 2(a-r) in Lugger et al.). 

NGC 1835 and NGC 1898 also appear to show power-law 
cusps and breaks, but at lower significance than NGC 2005 and 
NGC 2019. The profiles for these two clusters are shown in Fig. 
10(c) and (d) respectively, along with power-law fits to the central 
profile regions, and the best fitting EFF models. For each cluster, a 
power-law seems to model the central data better than an EFF pro- 
file, although for NGC 1898 the difference is small. The power law 
slopes are /3 = 0.45 for NGC 1835, and (3 = 0.30 for NGC 1898, 
and the break radii are 4'.' 2 ± 0'.'4 and 6'.' 3 ± l'.'2 respectively Pro- 
files for NGC 1835 have been published by both Elson & Freeman 
(1985) and Mateo (1987); however, while both show NGC 1835 to 
have a small core, neither find a power law region or a break. 

The three remaining clusters have incomplete profiles in their 
inner regions because of severe crowding. Even using short expo- 
sure photometry (Section 4.2.3) for NGC 1754 and NGC 1916, the 
incompleteness could not be overcome, while NGC 1786 had no 
supplementary images available. The profiles for these three clus- 
ters are shown in Fig. lO(e-g). Given the lack of data, we cannot 
assert even that these three clusters are PCC candidates; nonethe- 
less, they warrant further attention because they all appear to show 
breaks in their profiles at very similar radii to those observed for 
NGC 2005 and NGC 2019. Clusters with ordinary King-type pro- 
files should not show such breaks. Each of the three profiles is com- 
plete to several points further in than its break, and we attempted to 
fit power-law models to these points even though there are very few 
of them in each case - the idea being to explore the matter further 
by comparing potential power-law models to the best fitting EFF 
models. We obtained slopes of /? = 0.74 for NGC 1754, /? = 0.90 
for NGC 1786, and /3 = 1.17 for NGC 1916. Interestingly, these 
values are very close to those found for NGC 2005 and NGC 2019, 
especially given that the uncertainties for fits to only a few points 
(just three in the case of NGC 1754) are large. 

Profiles for both NGC 1754 and NGC 1786 are given by Ma- 
teo (1987); however he finds that neither show power-law regions 
or breaks. Similarly, Meylan & Djorgovski (1987) observe NGC 
1786, and they classify it as having a normal King-type profile. 
There is no published profile for NGC 1916. We find that for each 
of the three clusters, both the power-law profile and the EFF profile 
fit the data equally well around the break radius - hence we do not 
assert that these three clusters are PCC candidates. However, we 
believe that all three clusters warrant further high resolution obser- 
vations to obtain complete profiles. It is suggestive that not only do 
these clusters seem to show breaks in their profiles, but each ap- 
pears at least as compact as both NGC 2005 and NGC 2019, which 
are our leading PCC candidates, but which do not suffer as badly 
from crowding as do NGC 1754, NGC 1786 and NGC 1916. 

The results of the power-law fits for the seven clusters are 
summarized in Table 5, including the break radius for each clus- 
ter. On the basis of the above analysis, we conclude that NGC 2005 
and NGC 2019 are strong PCC candidates, while NGC 1754, NGC 
1786, and NGC 1916 deserve further detailed study. The status of 
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Table 5. Power-law slopes and break radii for seven potential PCC clusters. 



Cluster 


(") 


/3 


log Vbreak 
(") 


"^break 

(") 


NGC2005 
NGC2019 


3.63 ± 0.33 
3.61 ±0.24 


0.74 
0.71 


0.75 ± 0.02 
0.70 ±0.05 


5.6 ±0.3 
5.0 ±0.6 


NGC1835 
NGC1898 


4.76 ±0.26 
8.40 ±0.55 


0.45 
0.30 


0.62 ±0.04 
0.80 ±0.08 


4.2 ± 0.4 

6.3 ± 1.2 


NGC1754 
NGC1786 
NGC1916 


3.61 ±0.41 
5.45 ±0.38 
3.35 ±0.18 


0.74 
0.90 
1.17 


0.62 ±0.04 
0.75 ±0.04 
0.82 ±0.07 


4.2 ± 0.4 
5.6 ± 0.5 
6.6 ± 1.1 



Note: rc is the core radius of the best fitting EPF profile, from Table 4; jS is 
the slope of the best fitting power-law model; r^reak is the break radius 
observed from Fig. 10 with errors estimated by eye. 

NGC 1835 and NGC 1898 is unclear. While they both have central 
regions which are fit better by power-law models than EFF models, 
the power-law slopes are shallower than those tjTjically associated 
with PCC clusters. However, both show breaks at similar radii to 
NGC 2005 and NGC 2019. Physically, NGC 1835 is very similar 
in appearance to these two clusters, and we therefore consider it 
also a PCC candidate. NGC 1898 is a much more open cluster, and 
we speculate that it is perhaps on the verge of core-collapse. It is 
interesting that Olsen et al. (1998) find NGC 1898 to be ~ 3 Gyr 
younger than NGC 1835, NGC 2005 and NGC 2019. 

Given these results, we briefly consider the LMC old clus- 
ter population as a whole. In addition to the twelve clusters in the 
present sample (adding NGC 1466, NGC 1841, NGC 2210, NGC 
2257, and Hodge 1 1 to the seven discussed above) there are three 
more suspected old LMC clusters - Reticulum, which is definitely 
old and likely a member of the LMC based on its radial velocity 
(Suntzeff et al. 1992); and NGC 1928 and NGC 1939, which are 
definitely LMC members and likely to be old based on the spec- 
troscopic study of Dutra et al. (1999). Of the twelve measured in 
the present study, we suggest that 3 ± 1 are good PCC candidates. 
Of the three not measured in the present study. Reticulum is def- 
initely not (Suntzeff et al. (1992) estimate a core radius of ~ 13 
pc), and NGC 1928 and NGC 1939 are unknown as far as we are 
aware. Therefore we estimate that 20 ± 7 per cent of old LMC 
clusters are PCC clusters, a value which matches well the ~ 20 per 
cent estimated by Djorgovski & King (1986) for Galactic clusters. 
In addition, all the candidate LMC PCC clusters, as well as NGC 
1754, NGC 1786, NGC 1916, NGC 1928 and NGC 1939 (which 
all require further study) are in the bar region, at projected radii 
less than ~ 2.5° from the optical centre of the LMC. Similarly, 
Djorgovski & King find the Galactic PCC clusters to be centrally 
concentrated. 



5.1.4 R136 

R136 is the extremely compact, very young (t — 3 Myr) central 
cluster of the 30 Doradus H 11 region, which is the most luminous 
H II region in the local group. Because of its unique characteristics, 
we briefly describe it separately here. 

The surface brightness profile for R136 and the surrounding 
cluster NGC 2070 is clearly separable into two regions - an iimer 
region which is well fit by an EFF profile with rc = l'.'3 and 7 = 
2.43 to a radius r ~ 10" (see Fig. 6 and Table 4), and beyond this. 



a break to a shallower profile. This is consistent with the review 
of Meylan (1993) in which he finds a two-component King model 
necessary to fit his observed surface brightness profile. The first 
component in his model has core radius rc = l'.'3, equivalent with 
our measurement, and the second component a core radius = 
15". Following in this vein, we attempted a two-component fit to 
our profile. Fitting two EFF profiles, we found the best fitting core 
component to match our previously measured core profile, with the 
outer profile best fit by /i ~ 18.7, 7 ~ 2.6, and o ~ 39'.' 4, meaning 
a core radius rc ~ 33". This composite profile is plotted in Fig. 
1 1 (dashed line). It is clear that it provides a good fit to the entire 
profile, as opposed to our original one-component profile, which 
only fit the core. 

In the very outer part of our measured profile, we noticed a 
turn-down at low significance - this may be evidence for a tidal cut- 
off. If so, it would be more suitable to fit a King profile (Eq. 2) as 
the second component. Doing so, we found that once again the core 
was well fit by our original EFF model, with the King component 
having k ~ 18.2, rt ~ 130", and concentration c ~ 0.6, meaning 
a core radius rc ~ 33", which matches well the value from the 
two component EFF profile. The composite EFF-King profile is 
the solid line in Fig. 11. It is clear that it differs significantly from 
the two component EFF profile only in the very outer regions. 

Fitting two-component models to our overall profile did not al- 
ter our original core radius meanurement for R136, which suggests 
a degree of robustness to our result. There is considerable varia- 
tion in the previously published values of this core radius. Meylan 

(1993) finds rc = 1'.'3 from ground based data, while Elson et al. 
(1992) find rc = 0'.'5 from //5T observations. Malumuth & Heap 

(1994) obtain rc = C/.'QG from a mass-density profile based on HST 
data, and Hunter et al. (1995) find rc < 0"08, also using //5r ob- 
servations. Campbell et al. (1992) measure rc = 0'.'25 based on 
a pure power-law surface brightness fit, again derived from HST 
observations. Our measured core radius of rc = 1'.'3 is consistent 
with that of Meylan, but considerably larger than the other three. 
We reiterate however, that our measurement is an upper limit, due 
to the crowding in our images and also because our selected an- 
nulus widths mean that we are not equipped to measure core radii 
smaller than rc ~ l''5 accurately. Brandl et al. (1996) find con- 
siderable evidence for mass segregation in R136, and show that 
the derived core radius is therefore sensitive to the lower luminos- 
ity/mass cut-off for stars used in its calculation. For example, they 
find a core radius of rc = 0'.'48 using stars with m > I5M0, but 
increasing to rc = 0'.'97 using stars with m > 4Mq. This may go 
some way towards explaining the above variations. 

Similarly, there is some argument in the literature as to 
whether R136 might have undergone core collapse (see e.g., Camp- 
beU et al. (1992); Malumuth & Heap (1994)) - mostly the argu- 
ments are dependent on the core radius assumed for the cluster, 
and the relaxation time chosen to characterize the system. Camp- 
bell et al. calculate the relaxation time at the centre (their Eq. 2), 
and using their measured core radius (0.06 pc) and central density 
(5 X lO^Mopc"^) obtain an estimate of severalxlO"* yr - much 
shorter than the age of the cluster. They therefore suggest that R136 
is likely in an advanced state of dynamical evolution and has un- 
dergone core collapse. In contrast, Malumuth & Heap derive a core 
radius four times larger, and a correspondingly larger central relax- 
ation time - comparable to the age of the cluster. They also calcu- 
late the relaxation time at the median (half-mass) radius (their Eq. 
5) and obtain a similar estimate, therefore concluding that R136 is 
probably mass-segregated but not post core collapse (according to 
Biimey & Tremaine (1987) p527, core collapse occurs at 12 to 19 
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median relaxation times). Similarly, Brandl et al. (1996) show that 
dynamical mass segregation is likely to have occurred, at least for 
the most massive stars in the system, but that the median relaxation 
time is too long for core collapse to have taken place. 

Even though our resolution is not quite sufficient to provide 
more than an upper limit for the core radius of R136, our surface 
brightness profile can help shed some light on the issue. If R136 is a 
PCC cluster then its profile should show a cusp similar to those ob- 
served for the old clusters discussed in Section 5.1.3. Both Camp- 
bell et al. (1992) and Brandl et al. (1996) show that the profile of 
R136 can be fit by a power law (although in the case of Brandl et al., 
this is only for stars with m > 4QMq). We observe a small amovint 
of evidence for such a power-law cusp and a break. The central re- 
gion of the profile is plotted in Fig. 12, along with the best fitting 
power-law model and EFF model. We measure a slope /3 = 1.17 
with the break at r ~ 2'.' 6. This is only marginally consistent with 
the values of /3 measured for the old PCC candidates. From Fig. 12 
it is also clear that the EFF profile fits equally well, unlike for the 
best old PCC candidates. It is likely then that in the case of R136 
the apparent power-law structure in the central region of the pro- 
file is due to the brightest stars in the cluster residing in this space, 
consistent with the argvmient of Brandl et al. We can calculate the 
median relaxation time (in the form of Malumuth & Heap (1994), 
but see also e.g., Spitzer & Hart (1971); Spitzer (1987); Binney & 
Tremaine (1987); Meylan (1987)): 



trh — 



6.5 X 10'' 



ln(0.4iV) V m 



,.3/2 
h 



(Mtot\ 
V 105 ) 



1/2 



yr 



(7) 



where A'^ is the total number of stars and m is the typical stellar 
mass, Mtot is the total mass of the cluster and Vh is the median (or 
half mass) radius. Eq. 7 is essentially equivalent to Eq. 4.5 from 
Brandl et al. In the following Section (5.2) we derive estimates for 
the total masses of all the clusters studied in this paper. For R136 
we obtain Mtot ^ 2.5 x IO'^A/q, and using this and Eq. 11 we 
can solve for the half-mass radius, obtaining rh ~ 1.25 pc. Both 
this radius and the total mass are entirely consistent with the values 
obtained by Brandl et al. (3 x 10* M© and 1.1 pc respectively). 
Following Brandl et al., we adopt m ~ 0.5Mq and therefore esti- 
mate A' ~ 5 X IC*. Substituting these values into Eq. 7 provides us 
with trh ~ 2. 1 X 10* yr, very similar to the estimate of Brandl et al. 
{trh ~ 2.5 X 10* yr). For the most massive stars (i.e., m ~ 4OM0) 
we obtain trh ~ 2.1 x 10® yr, which is comparable to the age of 
the cluster. This suggests that for these stars, dynamical mass seg- 
regation has indeed had time to occur, but both calculated values of 
trh argue strongly against core collapse having happened. In con- 
firmation of this, we can estimate the central relaxation time (in the 
form of Campbell et al. (1992), but see also e.g., Spitzer (1987); 
Binney & Tremaine (1987); Meylan (1987)): 



:3xl0*fiM^oV 
\ m J 



0.06pc 



Po 



5 X lO^Mopc- 



1/2 



1.4 



ln(0.4Ar) 



(8) 



where Tc is the core radius and po the central mass density. In the 
following Section (5.2) we also estimate the central mass densities 
for all the clusters studied in this paper. Using our derived value 
of Po ~ 3 X lO^Afopc^^ for R136, and our upper limit for the 
core radius of Tc ~ 0.32 pc, we obtain t,.o ~ 1 x lO'^ yr for 
m ~ O.SMq, again considerably longer than the estimated age of 
the cluster. Given that tro oc Tc, the true core radius would have to 



be at least five times smaller than our upper limit in order to bring 
the central relaxation timescale down to the size required to make 
R136 dynamically old enough to be in the epoch of core collapse 
{tro ~ 1 X 10^ yr). 



5.2 Luminosity and mass estimates 

We can use our measured structural parameters to obtain luminosity 
and mass estimates for each cluster. Eq. 3 must first be deprojected, 
which is done by means of an Abel integral equation (EFF87; see 
also Biimey & Tremaine (1987), Section 4.2 and Appendix 1.B.4) 
for the luminosity density j (r ) of the cluster in question. We obtain: 



j{r) = jo 



(7+1) 



(9) 



which has the same functional form as the surface brightness pro- 
file n{r) but with index 7 + 1. In Eq. 9, jo represents the central 
luminosity density and is given by 



Jo 



Mor(^) 



(10) 



where F is a standard gamma function. To obtain the measured 
enclosed luminosity L as a function of radius, we integrate Eq. 9 
within a cylinder of radius r along the line of sight, since this is the 
relevant observational quantity. This gives: 



L(r) = 47rjo 



1 + 



+ t 



(t+1) 



d£ dx 



(11) 



where i is the radial variable and x the line-of-sight variable. Eval- 
uating the integral for r = rm, the maximum radial extent of the 
measured surface brightness profile, gives: 



Lm — 



( If '■ 

(a — a' {a 



2 , 2 N- 



(7-2) ' 



(12) 



By taking the limit Vm — > 00, we can also obtain an estimate for 
the asymptotic cluster luminosity Loo '■ 



2ivnoa 
7-2 



(13) 



provided 7 > 2, otherwise the limit is divergent. 

The calculated values for jo. Loo and Lm are listed in Table 
6. When calculating Lm we take from Table 4. In order to use 
Eq. 10, 12 and 13, we must have /xo, a and rm in physical units. As 
always, to convert to parsecs we take the LMC distance modulus 
to be DM = 18.5, which implies a scale of 4.116 arcsec pc~^. 
This takes care of a and r,n, but to convert /lo to Lq pc~^ is more 
complicated. We first need to know the V555 magnitude of the sun. 
We assume an absolute standard magnitude of V = +4.82 and a 
standard colour B — V = +0.65 and combine Eq. 7 and 8 from 
Holtzman et al. (1995b): 



WFPC2 = SMAG - Ti,Fs x SCOL 

—T2,FS X SCOL^ + Zfg — Zfs 



(14) 



where the notation is as in Holtzman et al. (1995b). Using our as- 
simied standard solar magnitude and colour, Zfg from Table 6 of 
Holtzman et al., and Ti,fs, 22, fs and Zfs from Table 7 of Holtz- 
man et al., we obtain U,>J^r, = +4.85. Alternatively, using Eq. 11 
and 12 from Dolphin (2000b) provides an identical result. The cen- 
tral surface brightness for each cluster must be corrected for absorp- 
tion, and we use E{B — V)= 0.10, which is a reasonable average 
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in the direction of the LMC (cf. EFF87). R136 however, has sig- 
nificantly higher reddening, with E{B — V)= 0.38 (Hunter et al. 
1995). Armed with these numbers, we use the relation 

log Mo = 0.4(Vb®5 - /i666(0) + DM 

+-i.lE{B - V)) + log(4.116^) L© pc"^ (15) 

to convert /io to physical units. 

We calculate mass and density estimates for each cluster by 
multiplying the appropriate luminosity equations by the average 
mass to light ratio MjLy for the cluster in question. The density 
p{r) then corresponds to Eq. 9, with the central density po obtained 
from Eq. 10. Similarly, the mass Mm inside Vm is given by Eq. 12 
multiplied M jLy and the asymptotic mass M^o is given by Eq. 
13 times this ratio. 

To estimate MjLy for each cluster, we use the evolution- 
ary synthesis code of Fioc & Rocca-Volmerange (1997) (PEGASE 
v2.0, 1999). This code determines the integrated properties of a 
synthetic stellar population as a function of time, using libraries 
of isochrones and stellar spectra to make the detailed calculations, 
and accounting for factors such as the IMF of the system, a range 
of metallicities, and any ongoing star-formation. We proceed using 
the simplest possible model - a population of stars formed simulta- 
neously in one initial burst and with the same metallicity - presum- 
ably a fairly good approximation to the formation of a rich stellar 
cluster. For the initial burst, we assume the IMF of Kroupa, Tout & 
Gilmore (1993) over the mass range 0.1 to 12OM0. There are four 
available initial abundances which cover the metallicity range of 
the cluster sample -Z = 0.0001 ([Fe/H] w -2.25); Z = 0.0004 
([Fe/H] w -1.65); Z = 0.004 ([Fe/H] w -0.64); and Z = 0.008 
([Fe/H] « — 0.33). The M/ Lv values for each metallicity are plot- 
ted as a function of cluster age in Fig. 13. Although the mass-to- 
light ratios obtained from these models are purely theoretical, they 
do agree well with observations (see e.g., Parmentier & Gilmore 
(2001)). From Fig. 13, it is clear that for most of the evolution, 
the M/Lv ratio is relatively insensitive to the chosen abundance. 
For completeness, we adopt the metallicity for each cluster based 
on the estimates in Table 2. Because of the insensitivity of the cal- 
culation to the selected abimdance, we are confident in using even 
those abundances which were calculated as averages rather than be- 
ing taken directly from the literature - it is sufficient to be able to 
differentiate between "metal poor" ([Fe/H] « —2.25) and "metal 
rich" ([Fe/H] w —0.33). We then use the age estimates from Table 
2 to obtain a mass to light ratio from the appropriate evolutionary 
synthesis model. The adopted abundances, and values for M/Lv, 
Po, Moo and Mm are listed in Table 6. 

The values of Loo and Moo are intended to provide reasonable 
estimates of total cluster Ivuninosities and masses in the absence of 
tidal limit information. There is very little difference in the result 
obtained from using Eq. 12 with Vm several hundred arcseconds (a 
reasonable tidal cut-off), and the result obtained from Eq. 13. We 
provide estimates of Lm and Mm to show the Ivuninosity and mass 
over the radial range we actually measure, and as such, these pro- 
vide reliable lower limits for total cluster luminosities and masses. 
For most clusters there is not a great difference between these val- 
ues and the asymptotic values; however, for clusters with 7 « 2, 
the value of Loo (and Moo) can become unreasonably large (e.g, 
NGC 1856), and the extrapolation r — > 00 may not be justified. 

The presence of saturated stars on some images, and the impo- 
sition of brightness limits on others (see Sections 4.2.3 and 4.2.4) 
causes some luminosity (and mass) to be unaccounted for in sev- 
eral surface brightness profiles (that is, the zero-point (1555 (0) is 



fainter than would be expected were the brightest stars included in 
the profiles). However, this effect is not particularly significant. It 
is most pronounced in the yoiuigest clusters, where bright stars are 
relatively massive. At worst we miss ~ 30 stars from any young 
cluster, which results in ^ GOOM© being neglected. Addition of 
this missing mass gives totals still within the quoted uncertainties 
for such clusters. The errors listed in Table 6 represent the uncer- 
tainties in the calculated parameters due to the random errors in the 
values of /io, 7 and a. We do not account for any systematic er- 
rors such as those introduced by any saturated stars or those which 
might be present in the calculated values for M/Ly . 



6 THE CORE RADIUS VS. AGE RELATIONSfflP 

With a large sample such as this, there is the opportunity for a full 
statistical analysis of the parameters presented in Tables 2, 4 and 6, 
in search of correlations and physical insight into the overall prop- 
erties of the LMC cluster system. Such an analysis is beyond the 
scope of the present paper and will be presented in the future. We 
do however observe one particularly noteworthy trend. 

When core radius is plotted against age for all clusters in the 
sample (Fig. 14), a clear relationship exists between these two pa- 
rameters - namely that the spread in core radius increases signif- 
icantly with increasing cluster age. Without exception, all of the 
young clusters have compact cores with rc ^ 2.5 pc, while for 
clusters older than log t ^ 9 the full range fc ~ — 8 pc is cov- 
ered. This trend has previously been discovered and discussed by 
Elson et al. (1989), E91 and Elson (1992), using a combination of 
literature data and ground-based measurements. 

The key question is whether this observed relationship is in- 
dicative of true structural evolution in LMC clusters as they grow 
older, or is merely the result of a secondary correlation, a product 
of the reduction process or a selection effect (i.e., we systemati- 
cally missed all the young, low surface brightness clusters from our 
sample). The fact that we confirm exactly the result of Elson, using 
a larger (and distinct) sample, space-based observations, and uni- 
formly selected and reduced data is argument against the relation- 
ship being a reduction artifact, or due to a selection effect. In addi- 
tion, we did not choose our young clusters in a systematic fashion, 
but rather simply used all those available in the WST archive - in ef- 
fect, a random sample. This adds weight to the argument against a 
selection effect being responsible for the observed upper envelope; 
however we cannot guarantee that there are not very low surface 
brightness young clusters with large cores in the LMC system. It is 
likely however, that such diffuse clusters would not remain bound 
for long, and so their presence (or otherwise) does not affect the 
discussion and conclusions below. 

Given that the radius-age relationship does not seem to be the 
product of data reduction or a selection effect, E91 discusses the 
possibility that a correlation between mass and core radius and be- 
tween mass and age might be responsible - that is, if both more 
massive clusters (with larger cores) and less massive clusters (with 
smaller cores) were formed in the past, but only less massive clus- 
ters (with smaller cores) formed recently, it might seem that core 
radius evolves with age. However, Elson observed no such correla- 
tions in her sample of 10 young clusters. 

Similarly, we plot Moo against age, and Moo against rc for 
our (much larger) sample, to look for correlations. These plots are 
Fig. 15 and Fig. 16 respectively. While Moo is not always a reli- 
able extrapolation to the true mass of a cluster, we use it here as a 
convenient parameter for eliminating the age bias in Lm and Loo 
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Table 6. Luminosity and mass estimates calculated using the structural parameters from the best fitting EPF profiles. 



Cluster log/i()° Adopted Adopted log jo log Loo log Lm logpo log Moo log Mm, 

(Lopc-2) [Fe/H] M/Lv (^oPC-^) (Lq) (Lq) {Mqpc-') (Mq) (Mq) 



NGC1466 


3. 


,06 




0, 


,02 


-2. 


,25 


2, 


,80 


2.33 ±0.05 


4.86 ±0.11 


4.80 ±0.09 


2.78 ±0.05 


5.31 ±0.11 


5.25 + 0.09 


NGC1651 


2. 


,54 


± 


0, 


,02 


-0. 


,33 


1. 
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4 70+0.10 
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NGC1711 


3. 


,38 
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.02 
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,64 
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NGC1718 
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" Corrected for reddening using E{B — V) = 0.10. 

^ Luminosity and mass calculated using the structural parameters derived for the inner regions only - that is, those parameters listed in Table 4. Given the 
break to a shallower profile at r ~ 15", we likely under-estimate Loo and Moo - the asymptotic limiinosity and mass including NGC 2070. 
Corrected for reddening using E{B — V) = 0.38 (see text). 
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(i.e., young clusters are intrinsically more luminous than old clus- 
ters), as well as any observational bias present in Mm (i.e., a cluster 
with 7 ~ 2 will have more mass at large radii, where we did not 
observe, than a cluster with 7 ~ 4, and consequently will likely 
have a correspondingly small Mm)- In Fig. 15, we see no signif- 
icant correlation of mass with age. It is true that there are a few 
young clusters with low masses, and that all the old clusters seem 
to have larger masses; however we attribute this to the fact that the 
low-mass young clusters will have faded and probably dispersed 
by age 10^" yr - explaining why there are no low-mass old clusters 
(a detailed discussion of this effect is presented in Meylan (1993)). 
There is no evidence for any significant intrinsic difference between 
the masses of the largest newly formed clusters and the old clusters 
in the sample. In Fig. 16 we also see no evidence for a strong cor- 
relation. While it is true that those clusters with large core radii all 
seem to be relatively massive, they are not more massive than most 
of the clusters with smaller core radii. 

We therefore conclude, like Elson, that the trend observed in 
Fig. 14 represents real evolution in the structure of clusters as they 
grow older. Given this, it is useful to return to Fig. 14 and exam- 
ine closely the distribution of clusters, which seems to exhibit a 
bifurcation at around several x 10* yr. This bifurcation grows into 
a large separation by 10 Gyr. It is possible that the apparent dearth 
of older clusters with core radii in the range 3 — 6 pc could be 
a small-sample effect. A simple Kohnogorov-Smimov test on the 
distribution of core radii for all clusters older than 10 Gyr, shows 
that they are not drawn from a uniform distribution, at about the 99 
per cent level. If we increase the sample size by including all clus- 
ters older than 1 Gyr, the significance increases to better than 99.5 
per cent. The bifurcation is therefore statistically significant. 

The lower branch in the radius-age diagram represents the 
standard picture of (isolated) cluster evolution, where a newly 
formed cluster suffers significant and rapid mass-loss due to stellar 
evolution, causing the cluster to expand. The rate of expansion is 
regulated by the initial mass spectrum present in the cluster (i.e., 
the slope of the mass function) (see e.g., Chernoff & Weinberg 
(1990)). The more heavily weighted the high mass end of a cluster's 
mass fimction is, the more severe the early mass loss. As the clus- 
ter grows older and evolves, the expansion slows and is eventually 
reversed by core collapse. E91 present Fokker-Planck models of 
evolving clusters with different IMF slopes, plotted over the equiv- 
alent of Fig. 14. Clusters with IMF slopes similar to that for the 
SalpeterIMF(a; = 1.35, where the IMF is given by $ oc m'^^^"^) 
and steeper follow tracks which match well the shape of the lower 
branch in Fig. 14. 

There is a noticeable scatter about this sequence, but its cause 
is not clear. Because the IMF slope governs the rate of expansion, 
IMF variations could be responsible. However, more and more ev- 
idence is pointing towards the universality of the IMF ~ Gilmore 
(2001) provides a detailed review, and more recently de Grijs et 
al. (2002c) have shown that a sample of six LMC clusters (widely 
scattered on Fig. 14) must have had very similar IMFs. It therefore 
seems likely that a combination of other factors is instead responsi- 
ble. Even the youngest clusters in the present sample have a spread 
in core radii - this probably reflects the different formation con- 
ditions of each cluster (E91), as well as the spread in initial mass 
(Fig. 15). Variations in initial mass would cause a spread in the 
lower sequence - it may be that this alone could reproduce the re- 
quired scatter. A second possibility is that the scatter is a result of 
the cluster sample being spread over a considerable range of dis- 
tances from the LMC centre (Table 2). Gilmore (2001) shows that 
tidal forces are at least as significant in the comparative evolution 



of Galactic globular clusters as IMF variations between clusters. 
The spread in the lower sequence of Fig. 14 might simply reflect 
the spread in orbital radii between clusters. We are currently em- 
ploying jV-body simulations to explore these two mechanisms in 
detail, to measure their relative effects, and to sec whether cither 
(or a combination of both) is able to reproduce the lower region of 
Fig. 14 for a suitable cluster sample. 

The fact that the models of isolated evolving clusters presented 
by E91 match the lower sequence of Fig. 14 so well adds plausi- 
bility to our argument that the radius-age plot tracks the physical 
evolution of LMC clusters. This renders the clusters which follow 
the upper branch especially intriguing. What makes these clusters 
different - why do they diverge from the standard sequence beyond 
a few X 10* yr, to finish in the upper right of the plot? Possibly, these 
clusters are on the verge of dissolution - after a certain time, mass 
loss can be so significant that rather than undergoing core collapse, 
a cluster simply dissipates. However, the upper sequence covers the 
timescale 10* to 10^'' yr, which is far too long for dissolution. In 
addition, Fig. 16 shows that clusters with large core radii are all 
relatively massive, whereas one would expect a cluster dissolving 
due to mass loss to be less massive. Upper sequence clusters are 
therefore probably not undergoing dissolution, but are expanding 
for another reason. Given this, wc suggest two ways in which up- 
per sequence clusters might be distinguished from the "standard" 
clusters - either they have (or have had) radically different stellar 
populations, or they have been influenced by some external process 
(i.e., they are not isolated). 

Examples falling into the first category include very flat IMF 
slopes, and very large binary star fractions. It seems vmlikely that 
either could cause the degree of core expansion observed. As dis- 
cussed previously, the IMF slope for a cluster regulates its early 
expansion, but the rapid expansion phase does not last the entire 
lifetime of the cluster (as would be required for the cluster to evolve 
along the upper sequence). In addition, as shown by E91, a cluster 
with a very flat IMF suffers such extreme early mass loss that it 
becomes unbound after only a very short time {-^ 10^ yrs). Our 
earlier argument with respect to the universality of the IMF also 
holds here. Similarly, although the binary star fraction in a clus- 
ter is important in the context of halting core collapse, binary stars 
alone cannot drive the large scale core expansion we have observed. 
To extract energy from a cluster's binary stars requires close stellar 
encounters, and an expanding core radius implies lower stellar den- 
sity and therefore fewer encoimters in a given period. It is imlikely 
therefore that significant cluster expansion can be driven by binary 
stars over the required timescale, even given an exceptional binary 
fraction. 

Two examples which fall into the second category are the 
merger of double clusters and the effects of strongly variable, per- 
turbative tidal fields. Again, from simple physical arguments, it is 
difficult to see how either mechanism could produce the scale of 
core expansion which is observed without destroying the clusters in 
question. Double cluster pairs are relatively conmion in the LMC 
(Bhatia et al. 1991; de Oliveira et al. 1998; de Oliveira et al. 2000), 
and Bhatia & Hatzidimitriou (1988) have shown that it is statis- 
tically unlikely that they are all coincidental alignments ~ hence 
some are likely gravitationally bound and interacting. An interact- 
ing pair will merge after a relatively brief period (Bhatia 1990), and 
iV-body studies by de Oliveira et al. (1998; 2000) have shown that 
a binary cluster merger can result in a single cluster with a struc- 
ture which is stable after ~ 200 Myr and distinct from the structure 
of the original clusters. However, it is unclear whether large scale 
core expansion is a result - small scale expansion implies many 
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mergers would be required, and this is statistically very unlikely. 
In the case of tidal forces, it is uncertain as to what effects a tidal 
field can have on a cluster core. Tidal forces increase as the cube of 
the distance from the centre of a cluster and therefore affect mainly 
the outer parts. A tidal force strong enough to significantly alter a 
cluster's core would probably destroy the cluster. However, as with 
the merger of double clusters, no models have been explored in 
the context of core expansion, and so neither mechanism can yet 
be quantitatively eliminated. We are currently carrying out detailed 
TV-body simulations (e.g., WiMnson et al., in prep.) to explore this 
further. 



7 SUMMARY AND CONCLUSIONS 

We have compiled a pseudo-snapshot data set of two-colour obser- 
vations from the HST archive for a sample of 53 rich LMC clusters 
spanning the full age range 10^ to 10^° yr. The emphasis has been 
on trying to make this compilation and the subsequent reduction 
process as homogeneous as possible without sacrificing data in- 
tegrity. We have also compiled literature estimates for the ages and 
metallicities of these clusters, again trying to maintain consistency 
as far as possible. From the HSr observations, we have constructed 
surface brightness profiles for the entire sample and obtained struc- 
tural parameters for each cluster, including the core radius and 
power-law slope at large radii. Using these parameters we have also 
estimated the total luminosity and mass for each cluster. These data, 
along with the surface brightness profiles, are available on-line at 
http://www. ast. cam. ac. uk/STELLARPOPS/LMCjclusters/. 

The surface brightness profiles show a rich amount of de- 
tail, with young clusters in particular exhibiting bumps, shoulders 
and dips in their profiles. We see evidence for double clusters in 
our sample, as well as post core-collapse clusters. The PCC can- 
didates are especially interesting - the two best examples show 
clear power-law profiles at small radii, with slopes (3 ~ 0.7 in the 
log fj, — log r plane. Our sample covers twelve of the fifteen definite 
old LMC globular clusters, and we are able to estimate that 20 ± 7 
per cent of these clusters are PCC objects, matching the 20 per cent 
estimated for the Galactic globular cluster system. We have also 
shown that R136 requires a two component fit to its profile, and 
that it is likely not yet in a PCC state. 

If core radius is plotted against age for the entire sample, we 
see that while all the young clusters have compact cores, the spread 
in core radius increases with age, with the oldest clusters cover- 
ing the full range of core radii measured. We have argued that this 
trend reflects real evolution in cluster structure with age. The dis- 
tribution of clusters on the plot suggests a bifurcation at several 
hundred Myr, with most clusters maintaining small cores consistent 
with standard isolated globular cluster evolution, but with several 
moving to the upper right of the diagram and evolving large diffuse 
cores. We suggest that these clusters must be different to the "stan- 
dard" clusters, either by having exceptional stellar populations, or 
by being subjected to an external influence. We are currently em- 
ploying Af-body simulations to explore several physical processes 
which fall into these categories. 
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FIGURE CAPTIONS 

Figure 1. Observation geometry and area correction points for cluster NGC 1841. The centre of the cluster falls on WFC3 at pixel coordinates (417, 401) 
and is marked with a small cross. Annuli of width 4" have been drawn to radius 100" and the random points falling in the field of view are plotted for every 
second annulus. 



Figure 2. Completeness functions for NGC 2213 (left) and NGC 2005 (right). The function for the PC is denoted by the solid line, that for WFC2 by the 

long-dashed line, that for WFC3 by the short-dashed line, and that for WFC4 by the dotted line. The functions have been integrated over position and colour 
to be a function of V555 magnitude only. Shown above the completeness functions are the PC F555W images of the core of each cluster, with exposure times 
of 120 s (NGC 2213) and 500 s (NGC 2005). The images are slightly cropped to be approximately 730 pixels (33'.'2) on a side. It is clear from these images 
that the core of NGC 2005 is extremely compact, on a scale less than ~ 160 pixels, while NGC 2213 is well resolved even in its very centre. 



Figure 3. An example of the inclusion of information from a short exposure to alleviate significant incompleteness due to crowding, in this case for NGC 
2005. Our long exposure image of the core of this cluster is shown in Fig. 2. (a) Siuface brightness profiles for all four annulus widths, calculated from 
the original exposure (solid circles) and short exposure (open squares). The dotted line shows the approximate radius r^jeu at which the profiles start to 
deviate, (h) Completeness functions for the short exposure (solid line) and original exposure (dashed line) within r^i^^,. The completeness functions have 
been integrated over colour and position within rdeu to be a function of V555 only. The vertical dotted line marks the brightness limit Vum imposed on the 
short exposure photometry to remove the offset between the two surface brightness profiles. 



Figure 4. Surface brightness profiles for an EI'F model (solid line) and empirical King model (dashed Une). The profiles both have core radii of rc = 10" 
and central surface brightnesses of fi(0) = 17. For the EFF profile we have chosen 7 = 2.5 which is approximately the median value measured by EFF87 

for young LMC clusters, and this implies a = 11'.'62 to obtain the correct rc For the King profile we have chosen a concentration c = 1.5, implying 
rt = 316" and k = 16.93. The core and tidal radii are marked. The vertical dotted line shows the typical maximum radius for surface brightness profiles 
constructed in the present paper. 



Figure 5. Backgroimd fit for NGC 2005. Models of the form of Eq. 5 were fit to the profiles with annulus widths 3" (filled triangles) and 4" (filled circles) 
beyond r = 20". The two estimated parameters were a ~ 4" and fj,o ~ 16.9. The solid line shows the average of the two best-fitting models, used to 

determine the background level Vtg = — 2.51og(^ = 21.83, marked with an arrow. Errors shown are those determined from the light distribution (Section 
4.2.5) rather than the Poisson errors, which are necessarily calculated after the background subtraction but before the final fit (see text). Note the apparent 
imder-estimation of the scatter at large radii. 



Figure 6. Background-subtracted F555W surface brightness profiles for each of the 53 clusters in the sample. The four different annulus widths are marked 
with different point types: I'.'S width are crosses, 2" width are open squares, 3" width are filled triangles, and 4" width are filled circles. Error bars marked 
with down-pointing arrows fall below the bottom of their plot. The solid lines show the best-fit EI'F profiles. For each cluster the core radius rc is indicated 
and the best-fit parameters listed. When converting to parsecs, we assiune an LMC distance modulus of 18.5. 



Figure 7. Measured values of 7 plotted against those determined by EFF87 from measurements extending to r ~ 250" for ten young clusters. Errors in the 
j/-direction represent the random errors determined as described in Section 4.2.5 and do not include any estimation of the systematic errors which may be 
introduced by our background fitting algorithm. The dashed line is plotted for reference and indicates equality between the two sets of measurements. 



Figure 8. Measured values of Tc plotted against those determined from the following studies: EFF87 (filled triangles); Mateo (1987) (filled squares); Elson 
(1991) (not seeing corrected) (filled circles). The cross marks the value measured for NGC 1856 by Elson (1992). Errors in the x-direction are those quoted 
by the various authors, and in the y-direction are those determined as described in Section 4.2.5. The dashed fine is plotted for reference and indicates 
equaUty between the measurements. Radii are in arcseconds to account for the different distance moduh adopted in each study. 



Figure 9. Images of the six clusters described in the text of Section 5.1.2. Top left: Our WFC2 image of the centre of NGC 1850 (lower left) and NGC 1850B 

(centre). The image is approximately 70" (730 WFC pixels) on a side. Top centre: NGC 2153. Top right: NGC 2213. Bottom row, left to right: NGC 2004, 
NGC 2155, and NGC 2159. Each of the latter five is an ij-band DSS2 image, approximately 195" on a side. 



Figure 10. Central surface brightness profiles and power-law fits for the seven PCC cluster candidates in the present sample. On the left are shown the four 
clusters with complete profiles; on the right those three with incomplete profiles, but which seem to exhibit breaks. The four different point styles represent 
the four annulus sets, as in Fig. 6. The best fitting power-law models are shown, and the radii of the breaks to the power-law regions are indicated by arrows, 
with errors indicated. The slopes /3 of the power-law models are listed. Note that because of the magnitude scale the slopes of the models as plotted are 2.5/3. 
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Figure 11. Two-component fits for the complete surface brightness profile of R136/NGC 2070. The four different point styles represent the four annulus 
sets, as in Fig. 6. We show two different two-component profiles. The solid line represents an EI'F profile for the core (R136) and a King profile for the outer 
region, while the dashed line represents an B¥¥ profile in the core and an ¥F¥ profile in the outer region. Best-fit parameters for each of the two models are 
as discussed in the text. Note that the two only differ significantly at large radii. 



Figure 12. Power-law fit to the inner core of R136. The four different point styles represent the four annulus sets, as in Fig. 6. The best-fit EFF and power-law 
profiles are shown as is the power-law slope and break radius (arrow), just as for the PCC candidate clusters in Fig. 10. Again, because of the magnitude 
scale, the slopes of the power-law model as plotted is 2.5/9. 



Figure 13. Mass to light ratios predicted from the single-burst stellar population models of Fioc & Rocca-Volmerange (1997) (PEGASE v2.0, 1999), using 
the IMF of Kroupa, Tout & Gilmore (1993) over the mass range 0.1 to 12OM0, and calculated for the four available abundances which cover the range 
spanned by the cluster sample: [Fe/H] = -2.25 (long dashes); [Fe/H] = -1.65 (soUd fine); [Fe/H] = -0.64 (short dashes); and [Fe/H] = -0.33 (dotted 
line). 



Figure 14. Core radius vs. age for all clusters in the sample. Ages are the Uterature estimates from Table 2, and the core radii as listed in Table 4. The points 
marking clusters with radii smaller than ~ 1 pc should be considered upper limits, as discussed in Section 5. 



Figure 15. Asymptotic mass vs. age for all clusters in the sample. Again, ages are the literature estimates from Table 2, while the masses are from Table 6. 



Figure 16. Asymptotic mass vs. core radius for all clusters in the sample. Core radii are from Table 4, and masses from Table 6. 
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